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SECTION  I 
INTRODUCTION 


This  report  deals  with  the  results  of  an  investigation  to  determine  the 
feasibility  of  using  molecular  seeds  in  high  performance  working  fluids  such  as 
hydrogen  to  absorb  energy  incident  from  a  distant  laser  source  and  utilize  this 
energy  to  raise  the  heat  content  of  the  working  fluid  which  provides  thrust  for 
a  rocket.  The  ultimate  objective  of  a  laser  thruster  program  is  to  generate  a 
high  temperature-high  pressure  gas  which  is  predominantly  hydrogen  in  order  to 
utilize  the  low  mass  amf  therefore  high  specific  impulse  of  a  high  temperature 
hydrogen  rocket  propellant.  The  seed  molecules  are  required  to  absorb  the 
laser  energy  and  thereby  increase  the  gas  temperature  to  the  5000  to  6000°K 
range.  The  study  is  thus  part  of  the  beamed  energy  concept,  a  highly  attractive 
approach  to  several  rocket  applications. 

In  this  study,  focus  was  placed  first  on  calculating  analytically  and 
then  measuring  experimentally  the  absorption  per  cm  of  several  molecular  seeds 
at  temperatures  up  to  6000°K.  The  analytical  study  considered  molecular  seed 
candidates  for  both  the  CO  (CO,  l^O  and  0D  seed  candidates)  and  CO2  (HjO, 

HDO ,  D20  and  NH-j  seed  candidates)  lasers,  and  the  results  of  that  study 
are  presented  in  the  present  report.'  Obtaining  temperatures  as  high  as  6000*K 
in  a  molecular  gas  such  as  hydrogen  by  conventional  methods  such  as  shock  tubes 
is  extremely  difficult,  so  a  novel  heat  source  was  used.  This  source  was  a  hot 
dense  plasma  created  in  the  gas  mixture  under  study  by  pulsed  laser  induced  gas 
breakdown  and  then  sustained  by  the  focused  output  of  7  kW  CW  CO2  laser  for  a 
time  interval  on  the  order  of  ten  seconds.  This  time  limit  was  set  by  cell  and 
window  heating  and  is  not  a  fundamental  limit  to  the  laser  sustained  plasma 
lifetime.  The  cell  containing ^he  mixture  under  study  had  been  placed  in  one 
arm  of  an  interferometer  raodifred  in  such  a  way  as  to  permit  simultaneous 
determination  of  the  phase  shift  and  power  loss  for  radiation  passing  through 
the  heated  region  of  the  cell  from  the  interferometer  light  source,  a  third 
CO2  laser  operating  at  a  wavelength  of  interest  to  this  study.  Processing 
the  phase  shift  and  power  loss  information  resulted  in  obtaining  both  the 
magnitude  and  the  temperature  dependence  of  the  absorption  per  cm  of  the 
mixture . 

The  results  of  the  analytical  investigation  are  briefly  summarized  for  92 
mixtures  of  the  candidate  in  at  11  atmospheres  pressure.  Of  the  candidate 
molecules  for  coupling  to  CO  laser  radiation,  the  CO  molecule  performed  the 

_  O 

best,  exhibiting  an  absorption  per  cm  of  10  at  temperatures  as  low  as 
300°K  and  a  chemical  stability  which  permits  this  magnitude  to  be  maintained 
to  temperatures  in  excess  of  6000°K.  Optical  saturation  of  CO  was  found  not  to 
be  a  problem  at  laser  intensities  of  interest  to  the  beamed  energy  concept. 

The  molecule  D2O  was  found  to  be  the  most  effective  coupler  to  CO2  laser 
radiation,  having  an  absorption  per  cm  in  excess  of  10~  in  the  temperature 
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interval  between  4(X)°K  and  4800°K.  The  molecule  H70  was  found  to  couple 
effectively  to  both  the  CO  and  COo  laser,  its  absorption  per  cm  10  over 
the  temperature  range  between  400°K  to  3800°K  in  the  former  case  and 
13U0°K  to  4800°K  in  the  latter.  The  ammonia  molecule  is  a  very  effective 
coupler  only  at  temperatures  below  about  1000°K  due  to  thennal  decomposition. 

The  experimental  investigation  was  carried  out  only  at  COj  laser  wavelengths 
although  a  high  performance  CO  probe  laser  was  constructed  for  measurements  at 
5  pm. 


In  carrying  out  the  experimental  investigation  at  CO2  laser  wavelengths,  two 
interesting  results  were  obtained:  First,  for  the  binary  mixtures  Hj/l^O  and 
H2/D2O  the  measured  magnitude  of  the  absorption  per  cm  was  significantly 
higher  than  calculated  analytically  while  exhibiting  a  temperature  dependence 
consistant  with  that  determined  in  the  analytical  study.  Second,  it  was  found 
that  in  the  absence  of  any  plasma  whatsoever,  the  degree  of  energy  absorption 
by  the  l^/lUO,  anc*  *^2^^3  mixtLjres  from  the  focused  high  power  CW 

CO2  laser  was  sufficient  to  heat  the  mixture  locally  to  temperatures  as  high 
as  390()°K,  producing  a  visible  discharge  of  light  from  the  high  power  laser 
focal  region.  The  first  of  these  results  indicates  that  l^O  and  D2O  are 
indeed  very  attractive  molecular  seeds  for  coupling  laser  radiation  to 
rocket  thruster  working  fluids,  being  even  more  attractive  than  was  indicated 
by  the  analytical  study.  The  second  result  indicates  the  possibility  that 
there  exists  a  set  of  (seed  concentration  and  compos i t ion )/( laser  incident 
intensity)  conditions  for  which  the  laser  powered  rocket  thruster  can  self 
activate,  that  is,  bring  itself  up  to  a  temperature  region  of  high  absorbtivity 
to  the  incident  laser  radiation  without  the  need  of  any  mixture  preheating  with 
electric  sparks  or  any  energy  adding  device  other  than  the  incident  laser 
rad i at  ion  itself. 

In  the  following  sections,  this  investigation  is  described  in  detail, 
beginning  in  Section  IT  with  a  description  of  the  analytical  results.  Section 
III  is  a  description  of  the  experimental  set  up,  the  method  of  data  gathering, 
and  ttie  process  of  data  reduction.  The  results  of  the  measurements  are  pre¬ 
sented  in  Section  IV,  and  conclusions  drawn  from  these  measurements  as  well  as 
recommendations  for  future  work  are  given  in  Section  V. 
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SECTION  II 

ANALYTICAL  RESULTS 


Introduction 


This  section  summarizes  the  results  of  the  analytical  study  to  determine  the 
optical  absorption  properties  of  several  molecules  at  two  wavelengths,  namely  5.0 
pm  and  10.6  pm,  in  the  electromagnetic  spectrum.  The  molecules  considered  were  CO, 

H2O,  HDO,  DgO,  NH^  and  0D.  For  any  molecule,  the  absorption  per  unit  length  is 
the  product  of  the  absorption  coefficient  of  the  molecule,  k,  and  the  molecules' 
concentration  u.  This  study  considered  mixtures  of  the  above  named  molecules  in 
hydrogen  at  temperatures  extending  from  300°K  to  6000°K,  and  therefore  resolved 
itself  into  considering  separately  the  temperature  variation  of  both  k,  due  to 
molecular  energy  level  positions,  partition  function  values,  induced  emission,  and 
collisional  pressure  broadening,  and  u  due  to  density  decrease  and  shifting  chemical 
equilibria  with  increasing  temperature. 

The  following  subsection  is  a  general  discussion  of  the  analytical  techniques 
used  to  study  the  temperature  dependence  of  both  k  and  u.  In  the  former  case  two 
general  approaches  are  used:  l)  a  line  by  line  summation  over  the  known  optical 
transitions  of  the  molecule  in  the  wavelength  region  of  interest  and  2)  an  empirical 
method  useful  for  those  situations  in  which  the  transitions  are  too  numerous  and  too 
complicated  to  permit  rapid  line-by-line  calculation  of  k.  The  former  method  is 
amenable  to  all  the  molecules  studied  at  temperatures  below  about  1500°K  as  both  CO 
and  0D  are  simple  enough  in  structure  to  permit  an  easy  all  inclusive  treatment  of 
the  transitions  of  interest  while  for  the  other  three  molecules,  tables  of  transi¬ 
tions  and/or  energy  levels  permit  calculation  of  k  so  long  as  the  temperature  is  not 
so  high  that  transitions  not  attainable  from  the  table  become  important.  The  line- 
by-line  calculation  is  relatively  useful  at  the  lower  temperatures  as  it  provides  a 
detailed  description  of  the  wavelength  dependence  of  k,  permitting  selection  of  an 
effective  laser  line  for  coupling  energy  into  the  molecule.  At  higher  temperatures, 
the  increase  in  spectral  line  density  due  to  excited  state  transitions  and  band 
overlap  causes  an  increase  in  the  accuracy  of  the  empirical  band  formulae  which 
provide  the  average  value  of  k,  denoted  k,  over  a  specified  wavelength  interval. 

The  results  of  the  application  of  these  techniques  to  the  above  mentioned 
molecules  are  described.  In  each  case  the  equilibrium  concentration  of  the  molecule 
is  calculated  and  combined  with  the  calculated  value  of  k.  Where  warranted,  a  dis¬ 
cussion  of  the  degree  to  which  equilibrium  is  established  in  the  mixture  under 
practical  conditions  is  included.  In  addition,  for  CO,  the  extent  of  optical  bleaching 
of  the  molecule  by  the  coupled  radi  on  field  is  calculated. 
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The  findings  of  the  study  indicate  CO  as  the  coupling  molecule  of  choice  at 
9  sm  wavelength  and  DgO  the  molecule  of  choice  at  10.6  urn. 

Analytical  Methods 

.Equilibrium  Chemistry 

For  each  of  the  candidate  molecules,  calculation  was  made  of  the  equilibrium 
chemical  composition  in  a  mixture  initially  containing  ten  atmospheres  of  hydrogen 
oil  one  atmosphere  of  the  candidate  molecule.  The  calculations  were  carried  out  for 
constant  total  pressure,  eleven  atmospheres,  at  temperatures  up  to  6000°K.  The 
calculations  were  made  using  a  computer  program  described  in  detail  elsewhere. ^ 

The  density  unit  used  throughout  this  report,  the  amagat ,  is  the  density  of  an 
ideal  gas  of  one  atmosphere  pressure  at  273°K,  2.69  x  lO1^  molecules  per  cubic  centi¬ 
meter. 

Absorption  Coefficient  Computation 
_Line^_by_-L_in_e_Cal£ul_at_ion 

In  this  method,  the  value  of  the  absorption  coefficient,  k,  at  frequency  cu 
is  expressed  as  a  sum  over  the  contributions  at  w  of  optical  transitions  centered 
at  frequencies 


k(u>)  =  I  kjlw.Wj) 


(1) 


where  the  frequency  of  the  ith  transition  and  k^(u>,(jjj)  is  the  contribution  of  this 
transition  to  k(a>).  The  value  of  k^w,^)  is  the  product  of  the  line  intensity,  Sit 
for  the  transition  and  its  associated  line  factor,  L^,  which  for  the  high  (10  atmos¬ 
pheres  hydrogen)  pressures  of  interest  to  this  study,  takes  the  form  peculiar  to 
collisions!  line  broadening: 


L,(oj,u>j,yj)  =  yj/(7r(yjZ  +  (wuii)2))  (2) 

In  this  Expression  y ^  is  the  line  width  of  the  transition  due  to  collisions  both 
v/ith  other  absorbing  molecules  and  with  other  molecules  in  the  mixture. 
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For  a  molecule  possessing  a  permanent  lipole  moment  M0  and  undergoing  a 
transition  from  a  rotational  angular  momentum  state  characterized  by  the  quantum 
number  J"  to  one  characterized  by  J'  the  value  of  S^,  denoted  by  S^(J",  J',  id)  for 
this  pure  rotational  transition  is  given  by  2 


S,(J,i,Ji,cj)=  s'?"  exp(-  E#/kT)(|-  exp(-hcuij/kT))  (3) 

3ncQlT)  J 


where  Eg,  is  the  energy  of  the  molecule  relative  to  J"  equaling  zero  in  the  lowest 
vibrational  state,  Q^T)  is  the  internal  energy  partition  function  of  the  molecule, 
and  SJ*,  is  the  line  strength  of  the  transition  and  depends  on  the  values  of  J"  and  J* 
Where  the  molecule  changes  its  vibrational  as  well  as  its  rotational  level,  is 
denoted  by  S^l”,  J’,  v"):3 


SBuh(v"+i)Sj* 

Si(j",J,,v",a>)=  ~jT-j - exp(-EJj/kT)(l-exp(-hcwj/kT))  (4) 


w0 


where  v"  is  the  lower  level  vibrational  quantum  number,  is  the  energy  separation 
between  the  levels  (v",  J"  =  0)  and  (v"  +  1,  J*  =  0)  and  Sg,  called  the  band  strength 
is  the  sum  of  S^CJ",  J’,  v",  w)  over  all  J"  and  v".  For  the  simple  molecules  CO  and 
OD  the  quantities  Ej,  and  S^.i  are  calculated  relatively  easily  and  the  use  of  the 

Line  oy  line  method  to  calculate  for  these  molecules  will  be  described.  Also, 

tabulations  of  wj_,  and  Sj»  exist  for  H20,  HDO,  D20  and  NH^  so  that  highly  resolved 
absorption  spectra  were  calculated  and  the  results  of  those  calculations  will  be 
given.  The  fractional  power  loss  per  unit  length  calculated  by  line-by-line  calcula¬ 
tions  is  denoted  by  a. 

Band_Mo_del  _C al cul_ati on s_ 


As  can  be  seen  from  the  above  discussion,  line  by  line  calculations  are  par¬ 
ticularly  useful  when  there  are  either  relatively  few  transitions  contributing  to 
the  value  of  k(w)  or,  if  this  is  not  the  case,  the  values  of  and  are  readily 
calculable  from  simple  formulas.  As  will  be  seen,  the  latter  proves  to  be  the  case 
for  simple  molecules  as  CO  and  OD  but  not  for  more  complex  ones  as  H2O,  HDO  and  D20 
and  NHo.  For  the  latter,  line  by  line  calculations  may  be  done  at  relatively  low 
temperatures  using  values  of  Et,  wi  and  obtained  from  experimental  data  and 
compiled  in  tabular  form  by  a  number  of  workers.  At  higher  temperatures  however, 
transitions,  originating  from  higher  lying  energy  levels  not  appearing  in  the 
tables  become  important  and  the  value  of  k(uj)  calculated  using  the  line  by  line 
method  is  too  small.  Since  for  these  molecules  the  necessary  additional  information 
needed  can  not  be  readily  obtained,  resort  is  made  to  band  models  in  yhich  the 
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measured  absorption  coefficient  centered  at  w  and  extending  over  a  given  frequency 
interval  to  either  side  of  w  is  set  equal  to  the  average  of  k(ou)  over  the  interval 
inside  which  the  line  intensity  and  the  interline  spacing  are  assumed  to  have  given 
mathematical  forms.  This  study  uses  the  work  of  Ludwig^*1*  in  which  the  spectral 
lines- of  H2O  in  both  wavelength  regions  of  interest  are  assumed  to  be  randomly  spaced, 
and  the  probability  of  a  line  having  intensity  is  assumed  to  have  the  form 

p(s,)=  s0_lexp(-s/s0)  (5) 

inside  the  spatial  interval,  25  cm~^in  Ludwig's  case,  for  which  measurements  were 
made.  For  this  model  the  experimentally  measured  value  of  transmittance,  T,  takes 
on  a  particularly  simple  form 

-in(T)  =  ku!(l  +  kul/4of^ 


where  k  is  the  value  of  ki  averaged  over  the  experimental  frequency  interval,  t  is 
the  experimental  pathlength  and  a  is  the  ratio  of  the  line  width  Yi  to  the  line  den¬ 
sity  d,the  number  of  absorption  lines  per  cm”l.  By  determining  T  at  a  given  tempera¬ 
ture  over  a  range  of  u,  one  can  determine  both  k  and  d,  assuming  a  knowledge  of  y^. 

From  his  measurements  on  HgO,  Ludwig1*  provides  values  of  k  and  d  for  the 
frequencies  of  interest  to  this  study  for  temperatures  up  to  3000°K.  As  the 
temperature  range  of  interest  in  the  present  study  extends  to  6000°K  and  the  species 
HDO  and  DgO  are  also  of  interest,  scaling  laws  were  developed  to  extend  Ludwig's 
results  to  higher  temperatures  and  the  other  isotopic  species  of  water,  and  these 
scaling  laws  are  now  described  in  general  terms  with  the  detailed  results  presented 
in  following  sections. 


The  scaling  laws  to  be  developed  are  derived  from  a  treatment  of  k  given  by 
Fenner. ^  in  this  treatment,  the  local  average  absorption  coefficient  k(J",  J',  v' , 
u> )  is  expressed  as  the  ratio  of  S(J”,  J',  v”,  w)  to  the  spacing  between  adajacent 
absorption  lines  which  for  a  simple  diatomic  molecule  is  simply  2B  where  B  is  the 
molecule's  rotational  energy  constant.  For  a  diatomic  molecule  in  a  electronic 
state  the  quantities  and  and  the  associated  values  of  o>  are  given  by 


and 


s  r-J+ 


ojj  =  cu0+  2B(j-fi) 


(7) 


Sj'1  =  J  (d) 

<jjj  -  oj0-  2BJ 
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For  a  pure  rotational  transition,  u)Q  is  replaced  by  zero  in  Fqs.  ^  and  8.  The 
partition  function  Q  is  factored  into  its  vibrational,  Qy,  and  rotational  Qr,  parts 

Ov=(i-exp(-hcw0/kT))''  ^9) 

0R=kT/hcB 


and  functions,  defined  as  k(J",  J',  v",  oj)  divided  by  (l-e""^Cll)/kT)  and  summed  over 
v",  can  be  written  in  the  following  form  for  rotational,  fr(tu),  and  rotation- 
vibration  ,  frv(w)  transitions  respectively 


(  8W3M02  \  /u>2  \ 

fR(w)  \^hc - J  V72/  exp(-<d2/<d’5 


do; 


where 


<d*  s(4kTB/hc)'/2 


Thus  it  is  seen  that  for  either  type  of  transition,  the  absorption  coefficient  data 
can  be  transformed  to  a  function  of  frequency  which  exhibits  a  maximum  at  oj  or 
equaling  w*,  which  is  seen  to  scale  as  (BT)1'  .  Furthermore,  in  the  case  of  fr(w), 
its  magnitude  at  w*  is  seen  to  be  independent  of  both  temperature  and  isotopic  species 
to  the  extent  that  MQ  is  species  independent.  In  contrast  g(i#)  is  expected  to 
increase  with  temperature  as  QV(T)  does,  but,  since  SB  itself  scales  with  i»0,  is 
dependent  on  isotopic  species  only  to  the  extent  that  QV(T)  is. 

In  the  following  subsections  it  is  assumed  that  {[(oj)  for  a  given  transition 
exhibits  the  same  general  behavior  as  described  above  for  k(J",  J',  v",  w)  so  that 
K(w)  may  be  scaled  as  described  for  temperature  and  isotopic  species  in  order  to 
obtain  ic  for  H2O,  HDO  and  D2  at  temperatures  up  to  6000°K. 
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Wavelength  Region:  5  urn 


Carbon  Monoxide 

£1  lemic  a  1  _C£a j"  u^  l  e  r  is  ti  cs 

The  equilibrium  composition  of  a  mixture  made  up  initially  of  ten 
atmospheres  Hg  and  one  atmosphere  CO  is  shown  as  a  function  of  temperature  in  Fig.  1. 
In  the  figure,  only  the  major  snecies  are  shown.  The  high  resistance  of  CO  to 
thermal  decomposition  at  high  temperature,  even  in  a  reducing  hydrogen  atmosphere, 
is  apparent  in  the  figure,  the  CO  fraction  having  been  reduced  by  only  50  percent 
at  6000°K.  At  the  lowest  temperatures,  less  than  1500°K,  the  equilibrium  situation 
is  apparently  dominated  by  the  water-gas  reaction  which  converts  CO  and  H2  to 
methane  and  water.  It  is  seen  that  in  this  temperature  range,  the  CO  fraction  is 
reduced  to  essentially  zero  if  chemical  equilibrium  is  attained.  If  it  is  assumed 
that  the  length  of  the  region  over  which  laser  radiation  is  to  be  absorbed  is  one 
meter  and  that  the  flow  rate  is  at  least  100  m/sec,  it  is  seen  that  the  time  avail¬ 
able  for  the  establishment  of  chemical  equilibrium  is  less  than  10  milliseconds, 
and  it  remains  to  determine  the  chemical  equilibrium  relaxation  time  for  the  system. 
It  is  known  that  reactions  involving  CO  and  H2  to  form  hydrocarbon  are  slow  and 
require  catalysts  at  temperatures  of  600°K,^  but  the  question  concerning  the 
degree  to  which  equilibrium  is  reached  at  temperatures  between  500°K  and  1500°K  in 
the  gas  phase  in  the  present  application  is  as  yet  unanswered. 

Spectroscopic^  £h ar ac te ri st £c £ 

Of  all  the  candidate  molecules,  carbon  monoxide  is  the  simplest,  with 
extensive  information  available  for  the  accurate  calculation  of  the  quantities  u)^, 
ujq  and  Et  in  Eq.  .  The  expression  of  for  S^m  is  given  by  Eqs.  7  and  ti  and  the 
value  of  Sg  has  been  determined  experimentally. ^  The  value  of  Q(T)  is  readily 
obtained  from  a  simple  formula®  incorporating  the  spectroscopic  constants  contained 
in  Ref.  9.  The  value  of  y  has  been  determined  for  CO  self-broadening  as  a  function 
of  J",10  and  the  value  for  H2  broadening  was  obtained  simply  by  multiplying  the 
elf-broadening  coefficient  by  the  square  root  of  the  ratio  of  the  CO-CO  reduced 
mass  to  the  Hp-CO  reduced  mass  under  the  assumption  that  the  difference  between  the 
two  coefficients  is  due  mainly  to  the  different  relative  velocities  for  the  two 

'Tiding  molecular  pairs.  Consistent  with  the  latter,  the  temperature  scaling  of  Y 
was  taken  to  be  T-1^  since  at  constant  pressure,  the  total  density  decreased  as 
T"  while  the  relative  velocity  of  collision  increased  as  T-^/^.  jn  calculating 
k(to),  the  summation  included  both  the  I'(J"-»J"-l)  and  R(j"-*-J"+l)  transitions  arising 
from  the  lowest  Uo  rotational  levels  in  the  first  ten  vibrational  levels  of  the 
molecule.  The  results,  are  presented  as  a  function  of  u>  and  T  in  Fig.  2  neglecting 
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chemistry  effects.  At  a  frequency  of  2012.73  cm”1,  the  value  of  a  at  300°K  is  seen 
to  be  dominated  by  contributions  from  the  wings  of  transitions  not  centered  at  this 
frequency,  for  the  lowest  energy  level  from  which  a  transition  occurs  near  2012.73 
cm”  ,  the  PQ(30)  at  2013. 31*  cm"1,  lies  1782.65  cm”1  above  the  lowest  CO  energy 
level  and  is  thus  very  sparsely  populated.  (The  notation  Po(30)  denotes  the  transi¬ 
tion  between  the  J  equals  30  level  of  the  lowest,  v  equals  zero,  vibrational  state 
to  the  J  equals  29,  a  P  transition,  in  the  first  excited,  v  equals  one,  vibrational 
state.)  As  the  temperature  increases  so  does  the  intensity  of  this  line,  particularly 
with  respect  to  its  neighboring  wing  background,  and  at  1000°K  the  value  of  a 
becomes  significantly  augmented  by  the  transition  which  appears  first  as  a 

weak  and  then  as  a  prominant  shoulder  on  the  Po(30)  transition.  In  addition  all  the 
peaks  appear  sharper  since  each  ^  is  decreasing  with  increasing  temperature.  Still 
other  contributing  transitions  become  important  as  the  temperature  is  increased  to 
6000°K,  the  P2(19),  the  P3(l3)  the  Pl*(7)  the  Rg(7),  RT(15),  the  Rq(25)  and  the  RQ(36) 
being  the  largest  single  contributors  in  their  bands  with  the  total  R  contribution 
being  about  an  order  of  magnitude  smaller  than  the  P  contribution.  In  Fig.  3  the 
temperature  dependence  of  k(2012),  in  units  of  cm”1  amagat”1  is  shown,  and  the 
absorption  per  cm  reflecting  the  temperature  variation  of  the  equilibrium 
chemistry,  is  shown  in  Fig.  U  along  with  a  curve  denoting  the  temperature  dependence 
expected  neglecting  chemistry  effects.  Included  in  the  curve  including  chemistry 
effects  is  a  contribution  to  the  absorption  due  to  the  presence  of  water,  to  be 
discussed  later.  Clearly,  HgO  is  the  dominant  absorber  in  equilibrum 
mixtures  below  1000°K.  However,  assuming  that  the  equilibrium  situation  is  not 
attained  until  significant  hydrogen  dissociation  is  achieved  at  temperatures  in 
excess  of  1500°K,  the  temperature  variation  of  CO  is  expected  to  follow  the  upper 
curve  at  temperatures  below  1500°K  and  the  lower  curve  thereafter  with  the  resulting 
absorption  coefficient  nowhere  falling  below  10”2  cm”1  in  magnitude. 

Optical  Saturat_ion_qf  C0_ 

The  essential  principle  of  rocket  propulsion  by  beamed  electromagnetic 
energy  is  that  the  energy,  once  absorbed  by  the  molecule,  be  rapidly  transferred 
from  internal  excitation  to  random  translational  motion  of  the  molecules  in  the 
medium  of  which  the  coupling  molecule  is  part,  with  this  transformation  being  accom¬ 
plished  by  collisional  relaxation  of  the  excited  molecule.  Specifically,  if  Ny  denotes 
the  number  density  of  molecules  in  the  v^*1  vibrational  level  of  the  molecule,  the 
change  in  Ny  with  time  will  be  determined  by  the  relative  rates  of  optical  and  col¬ 
lisional  excitation  and  de-excitation  of  the  molecule.  Mathematically,  this  may 
be  expressed  as  follows: 

Nv  -  N v-f| 

C 

(13) 

ri,VtV-t  {Nv-exp(-AEv%v_t/kT)Nv_,} 

+Z  rj<v  +t  v  |Nv+i  “  exp(— AEvv_|/kT)  nv| 
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where  the  first  and  second  terms  represent  transitions  between  the  vtfl  and  (v-l)bh 
and  (v+l)111  levels  respectively  caused  by  the  presence  of  a  radiation  field  of 
frequency u)  and  intensity  I  interacting  with  molecular  absorption  cross  sections 
°v-l,v  and  ov>v+l«  For  a  P  branch  transition11 

C  =  exp(-2JBhc/kT) 


The  second  pair  of  terras  in  Eq.  13  represents  transitions  caused  by  collisions 
with  other  molecules.  In  the  above,  T  represents  the  gas  kinetic  temperature,  while 
for  the  purpose  of  this  calculation,  an  effective  vibrational  temperature  Tv  is 


defined  such  that 


Ny 

Nv.| 


=  A-^expe-AE^/kTy) 


where  AEv+qjV  is  the  energy  difference  between  the  two  vibrational  levels,  neglecting 
rotation.  Invoking  the  harmonic  oscillator  approximation, 

^v.v+i1  (v+Oo'ai 

ffV-I.V  = v  *0.1 

Landau-Teller  selection  rules 

kV,V-l  =  Vkl,0 
kV+l,V  =(V  +  |)kl,0 

and  making  the  substitution 

B2exp(AEv+lv/kT), 


F.q.  13  can  be  written  to  express  the  steady  state  in  the  form 

Q  -  1  a0.'  1<70.l  Uv+i)-  (v  +  l)| 

hccu  l  hew  j  AC  I 


(iM 


-r. 


'10 


(v+  l) 
A 


V+J_l 

B  f 


and  the  problem  becomes  one  of  solving  this  equation  for  A,  thereby  obtaining  Tv 
and  the  sensitivity  of  the  latter  upon  I.  In  the  limit  of  zero  intensity  Tv  and 
T  are  seen  to  be  identical  while  in  the  limit  of  infinite  intensity,  Tv/T  is  seen 
to  equal  w/(2JB)  a  quantity  greater  than  unity,  reflecting  the  fact  that  without 
the  thermalizing  effect  of  collisions,  the  molecule  becomes  vibrationally  hot  in 
which  case  the  rate  of  stimulated  emission  tends  to  equal  that  of  absorption,  and 
the  molecule  thus  ceases  to  couple  to  the  radiation  field.  At  intermediate  values 
of  I  tht  degree  cf  decoupling  is  less  and  the  saturation  intensity,  Is  is  defined 
as  the  value  of  I  for  which  the  degree  of  optical  coupling  between  the  vth  and 
(v-l)th  levels  is  half  its  value  at  zero  intensity.  That  is,  Is  is  the  value  of 
I  at  which  A-  l/C  =  (B- 1/0/2 

The  value  of  r1Q  used  in  the  calculation  was  that  reported  by  Millikan  and  White12 
for  ten  atmospheres  of  H2  relaxant.  The  value  of  aol  is  the  value  of  k(w),  in  units 


19 


' 


R79-922895-25 


of  cm^  and  calculated  for  CO  using  Eqs.  U  and  2  for  L^(w^,  ui^,  y^,  divided  by 
the  stimulated  emission  factor  ( 1-exp ( -hcui^/kT) ) , 


°oi 


Si(j1;j-.oicui) 
w  Xj  (i-expf-hcojj/kr)) 


(15) 


and  aQ2  is  presented  in  Fig.  5  as  a  function  of  and  T.  It  is  seen  that  over 
the  entire  temperature  range  of  interest  to  this  study  the  value  of  the  maximum 
value  of  oQl  is  essentially  10  cm  ,  and  this  value  was  used  in  the  calculation, 
making  the  calculated  value  of  1=  a  lower  bound  on  the  true  value  in  the  case 
where  the  true  value  of  is  smaller  than  10“  7  cm  . 


Considerable  physical  insight  is  gained  by  making  several  simplifying  approxi¬ 
mations  with  the  resulting  value  of  Is  being  essentailly  identical  to  that  obtained 
by  solution  of  Eq.  lU  by  computer.  These  approximations  are  first  the  setting  of 
C  equal  to  unity  and  second  the  removal  of  the  v  dependence  of  each  term  in  Eq.  lh . 
With  this  done  it  is  found  that  Is  is  the  laser  intensity  for  which  the  rate  of 
excitation  of  the  first  vibrational  level  by  the  radiation  field  equals  the  rate 
of  deexcitation  of  that  level  by  collisions: 


Is  =  rl0  hcw/o-0l 


(16) 


In  Fig.  6  the  temperature  dependence  of  Is  for  a  mixture  of  one  atmosphere  CO  in 
ten  atmospheres  H2  is  shown  for  the  oQl  value  discussed  ab^ve.  It  is  seen  that 
for  the  gas  mixture  considered,  Is  greatly  exceeds  10  W/cra  ,  an  intensity  of  inter¬ 
est  to  the  beamed  energy  application,  over  the  entire  temperature  range.  It  is 
noted  that  the  large  values  of  Ig  obtained  are  due  entirely  to  the  effectiveness 
of  H2  as  a  CO  vibrational  quencher.  Were  He  is  to  be  used  in  place  of  H2,  the 
value  of  1^  can  be  expected  to  be  reduced  by  nearly  two  orders  of  magnitude  at 
300°K  and  approximately  a  factor  of  2  at  6000°K. 

h2o 

Ch  em  ic  al_C  haract  e.ri  s  ti  c_s_ 

The  chemical  composition  of  a  preequilibrium  mixture  of  one  atmosphere 
H20  and  ten  atmospheres  H2  is  shown  in  Fig.  7  as  a  function  of  temperature.  It  is 
seen  that  the  HgO  mixture  retains  its  initial  composition  until  a  temperature  of 
about  k000°  is  attained  at  which  point  dissociation  of  the  molecule  occurs. 
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Spe^ctro_sco£ic_  Characteristics 


Although  being  chemically  quite  simple,  the  H2O  molecule  is  extremely 
complex  spectroscopically.  Because  of  its  physical  structure  and  the  relative 
sizes  of  its  moments  of  inertia,  HpO  is  an  asymmetric  rotor  for  which  in  principle 

T»»  TO 

the  values  of  Ej ,  W3  and  Sj,  can  be  calculated  with  considerable  computational  effort. 
Moreover,  because  of  the  small  mass  of  the  hydrogen  atoms  in  the  molecule,  pertura- 
bation  of  the  spectroscopic  properties  of  the  molecule,  both  013  and  Sj,,  is  con¬ 
siderable,^  and  for  accurate  values  of  these  parameters,  resort  must  be  made  to 
tabular  compilations  where  such  exist. 


The  transitions  of  importance  at  5.0  ym  involving  HgQ  are  the  R  branch  of  the 
v p  as  well  as  both  the  P  and  R  branches  of  the  vp-vg  and  V3-V2  transitions.  The 
rotational  energy  levels  of  the  following  vibrational  levels  were  used  in  the  line 
by  line  calculation:  (000),  (010),  (020),  (100),  and  (001).  Each  rotational  level, 
designated  by  the -quantum  number  J,  has  2J+1  sublevels  each  designated  by  the  quan¬ 
tum  number  T,  and  the  extent  of  the  data  base  used  in  the  line  by  line  calculation 
is  indicated  by  the  magnitude  of  the  largest  J  included  in  the  energy  level  compil¬ 
ation.  For  the  (000)  level  the  rotational  sublevel  compilation  was  complete,  that 
is  all  2J+1  sublevels  were  available  in  tabulations  for  each  J,  for  J  values  less 
than  16  with  partially  complete  tabulations  available  up  to  a  J  value  of  22. 1 
For  (OlO),  the  tabulations  were  complete  up  to  J  equals  13  and  partially  so  up  to 
J  equals  20. ^  The  (020)  was  complete  up  to  J  equals  6  and  partially  so  up  to 
J  equals  11. ^  The  (100)  was  complete  up  to  J  equals  U  and  partially  so  up  to  J 
equals  8. -*-7  And  the  (001)  was  complete  up  to  J  equals  6  and  partially  so  up  to  J 
equals  10.^ 

Each  of  the  2J+1  sublevels  bears  in  addition  a  symmetry  designation  (++),  ( — ), 

(+-)  or  (-+).13  Because  of  symmetry  and  nuclear  spin  effects  the  value  of  S3  for 
transitions  originating  in  a  (++)  or  a  ( — )  level  includes  a  numerical  factor  of 
lA  while  for  a  transition  originating  in  a  (+-)  or  a  (-+)  level  includes  a  factor 
of  3A.  The  vg  and  transitions  are  both  known  as  type  b  transitions  so  that 

only  transitions  of  the  type  (•*■+)  -  ( — )  and  (+-)  -  (-+)  are  allowed.^  The  v  _ 

\>2  transition  is  a  type  a  transition  and  only  transitions  of  the  type  ( — )  _  (♦-) 
and  (■*■+)  -(-+)  are  allowed.1^  Finally  in  any  transition  the  absolute  value  of  J 
must  either  be  unchanged  or  change  by  no  more  than  unity. 


Given  the  available  known  energy  levels  and  the  selection  rules  described  above, 
all  the  allowed  transitions  between  these  levels  were  determined,  and  those 
occurring  between  energy  levels  separated  by  more  than  1965  cm--*-  and  less  than 
2065  cm-1  were  included  in  the  present  line-by-line  calculation.  As  indicated  by 
Eq.  it  is  first  necessary  for  each  transition  to  calculate  S3.  For  many 

of  the  stronger,  S3  greater  than  10-2'’  cm/molecule,  transitions  at  296°K,  this  cal¬ 
culation  was  already  carried  out  by  McClatchey,  et.  al,^  and  for  those  transitions 
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S^(T)  vas  calculated  from  the  expression 


S,  (296)0(296) 


S,W  = 


where  the  following  expression-1 

10 


e£(T-296) 
296  T 


I-  expUiCujj/kT) 
i  —  exp  (he  oj.  )/296 k) 


Q(l):l/2  I  exp  (VBC"/4kT)  Vrrk3T3/ABC  exp(-E/kT) 

YVV° 

is  used  to  evaluate  the  rotation  vibration  partition  function,  the  quantities  A, 

B,  C  and  E  depending  on  v^,  Vg  and  v^  according  to  the  simple  formulae  given  by 
Benedict  et  al.?^  Limiting  the  summation  to  10  for  v-p  Vg  and  v^  gave  convergence 
over  the  entire  temperature  range  of  interest. 

Many  of  the  possible  transitions  were  not  strong  enough  to  be  included  in 
McClatchey's  compilation  because  they  originate  at  levels  having  a  large  value  of 
Ej ,  causing  e-E&/k.T  to  be  small.  Such  transitions  may  become  important  and  even, 
if  their  line  strength  is  large,  dominant  at  higher  temperature,  and  the  inclusion 
of  such  lines  in  the  present  compilation  is  thus  desirable.  For  these  lines  the 
value  of  the  line  strength  was  calculated  using  Wacker  and  Pratto's  tables.^®  For 
the  R  branch  of  the  \>g  transition,  such  line  strengths  are  upper  bounds  because  of 
centrifugal  stretching  effects. Trial  calculations  using  known  strengths?  as 
comparisons  revealed  that  the  actual  strength  may  be  as  more  than  50  percent  lower 
than  calculated.  However  no  atterpt  was  made  to  correct  the  calculated  line 
strengths  for  centrifugal  distortion  effects.  Inclusion  in  the  present  compilation 
required  that  the  calculated  line  strength  exceed  10-4.  The  largest  possible  value 
of  Sj„  is  J+l.  With  the  line  strength,  E^ ,  ,  Sg  and  Q(296)  in  hand,  S(296)  was 

calculated  to  be  used  as  shown  in  Eq.  IT  in  calculations  at  higher  temperatures. 

The  values  given  by  McClatchey  et  al.?  for  y ^  pertain  to  line  broadening  by  Ng 
and  0g  in  the  atmosphere  and  are  note  relevant  to  the  present  study.  Instead,  recourse 
was  made  to  the  work  at  29o°K  of  Benedict  and  Kaplan^ ’^3, 2^4  which  gives  the  value  of 
Y-j  in  terms  of  J"  and  i",  the  quantum  numbers  defining  HpO  rotational  energy  levels. 

As  mentioned  previously  the  H?0-U?  values  of  were  converted  to  HpO-Hp  values  by 
i  ultiplying  by  /yNp-HpO/ijHp-HpO .  For  transitions  involving  (J",  t")  values  outside 
the  range  studied  by  Benedict  and  Kaplan,  the  value  of  y^  was  obtained  by  extrapola- 
ing  the  appropriate  curve  of  y^  vs  J"  at  constant  (J"  ♦  t")  to  the  desired  value 
of  J".  For  HpO  self  broadening,  y^  was  assumed  to  scale  as  T“?  and  for  Hp  broaden¬ 
ing  the  scaling  was  assumed  to  be  T“0*?.  The  results  of  the  line  by  line  calcula¬ 
tion  are  shown  in  Fig.  8  for  a  10:1  HgO  mixture  at  11  atmospheres  total  pressure 
in  the  temperature  range  300  to  6000°K.  In  the  frequency  internal  studied  the 
dominant  absorption  lines  arise  from  the  R  branch  of  the  band  due  to  its  large 
Gg  value  and  the  fact  that  all  its  transitions  arise  from  relatively  low  lying  energy 
levels,  the  value  of  a  decreasing  with  increasing  .  As  the  temperature  is  increased  to 
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about  2000°K  additional  transitions  become  important  and  more  fine  structure 
appears.  Apparently  the  spectrum  calculated  from  the  energy  level  data  base  is 
fully  developed  by  about  2000°K,  that  is, all  exp(-E  1/ kT)  are  on  the  order  of  unity, 
after  which  point  the  energy  level  compilation  is  inadequate  to  calculate  a,  and 
band  model  theory  must  be  used.  Used  as  a  tool  to  determine  desirable  laser  lines 
for  beamed  energy  coupling.  Fig.  8  reveals  that  the  region  around  1980  cm--*  to  be- 
relatively  transparent  at  300°K,  a  fact  consistant  with  McClatchey  and  Selby's2^ 
results,  while  becoming,  relatively  opaque  at  temperatures  exceeding  100CTK. 

As  outlined,  Ludwig's  data  was  reduced  to  give  the  function  g(w)  and  in  turn 
w*  and  a  factor,  referred  to  here  as  C  which  is  given  in  Eq.  11  by  the  quantity 
Qv(T)/u>o.  In  carrying  out  this  analysis  w0  was  set  equal  to  1600  cm”-'-  to 
facilitate  computation.  Figures  9  and  10  present  the  calculated  value  of  g(w) 
vs  u>  at  600°K  and  3000°K  respectively.  Also  included  are  curves  in  which  the 
parameter  set  (C,  w*)  has  the  values  (.111,  160)  and  (.163,  300)  respectively  for 
the  two  temperatures.  It  is  seen  that  the  single  expression  for  g(w)  given  in 
Eq.  11  is  a  good  representation  of  data  so  long  as  u>-l600  is  less  than  u>*.  At 
3000°K,  the  frequency  range  of  interest  to  the  present  study  does  lie  within  this 
range.  In  Fig.  11  the  parameters  C  and  w*  derived  from  fits  of  Ludwig's  data,  are 
shown  as  functions  of  T.  It  is  seen  that  w*  scales  essentially  as  which  is 

represented  by  the  drawn  straight  line,  and  that,  with  the  exception  of  the  point 
at  300°K,  C  is  given  by  0.098  ( 1-exp (-hcw0/kT) )-1,  where  the  cofactor  of  O.O98  is 
an  approximation  to  Q^T)  with  equaling  1600  cm-1.  Finally  it  is  noted  that^the 
value  of  0.098  corresponds  for  the  indicated  value  of  w0,  to  a  value  of  160  cm ~<L 
amagat”'-  for  S3,  compared  to  an  experimental  value  of  300  +60  cm"2  amagat"'-. 
Considering  the  simplicity  of  the  treatment  leading  to  Eq.  11  and  the  fact  that 
the  line  intensities  for  R  branch  HpO  transitions  are  expected  to  be  smaller  than 
calculated  due  to  centrifugal  distortions  effects,  this  value  of  Sg  is  in  not  too 
bad  agreement  with  experiment.  More  to  the  point  is  the  fact  that  the  analytical 
approach  used  has  provided  the  means  for  obtaining  the  quantity  E  in  Eq.  6  at 
temperatures  in  excess  of  Ludwig's  experimental  range  for  H2O  in  the  region  1965- 
2065  cm-1,  and  k  is  given  as  a  function  of  temperature  for  2000  cm-1  frequency  in 

Fig.  12  along,  for  comparison,  with  the  line-by-line  results  at  2012.5  cm”-*-  and 

2010  cm--*-  averaged  over  10  cm”-'-  to  either  side. 

From  Eq.  6  it  is  clear  that  the  value  of  a  .  n  cm  ^  is  Ru/(l  +  Ku/‘<a) 

where  u  is  the  HpO  pressure  in  amagats  and  i  is  one  centimeter.  Malkmus  provides 

the  value  of  d-1  as  a  function  of  T  in  the  frequency  range  of  interest  as  well  as 
the  line  width  coefficients  y ^  necessary  to  calculate  the  quantity  a  for  use  in 
the  above  expression.  It  is  found  that  for  the  gas  mixture  of  interest  (one  atmos¬ 
phere  U^.0,  ten  atmospheres  112),  Y  decreases  from  1.3  cm"-'-  to  0.21  cm-'-  as  T 
increases  from  300°K  to  6000°K  while  d-1  increases  from  .16  cm  to  bbOO  cm  over 
this  temperature  range  at  2000  cm-'-  so  that  (l  +  ku/^a)'^2  is  essentially  unity 
throughout  the  temperature  range,  and  acorresponds  to  ku. 
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FIG.  10 


ABSORPTION  COEFFICIENT  FUNCTION  g(W)  VSW  FOR  H20 

T  =  3000°K 


•  LUDWIG’S  DATA 


CM."2  AMAGAT- 


AMAGAT - 


R79-922896-25 


FIG.  12 
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In  Fig.  13  ku  is  given  as  a  function  of  temperature  at  2000  cm-1  for  both 
the  case  in  which  no  chemical  decomposition  of  HgO  is  assumed  to  occur  and  the  case 
in  which  the  H20  pressure  is  assumed  to  be  the  equilibrium  value  shown  in  Fig.  7. 
Examining  Figs.  12,  13  and  8  and  recalling  that  k  is  derived  from  measurements  which 
are  essentially  an  average  over  a  25  cm'1  interval,  the  results  of  the  two  cal¬ 
culations  seem  consistent  at  temperatures  below  2000°K  above  which  the  data  base  of 
the  line  by  line  calculation  becomes  inadequate  for  accurate  calculation  of  a.  It 
is  seen  that  the  absorption  per  unit  length  for  H2  is  on  the  order  of  10"2  cm-1,  for 
a  desirable  CO  laser  line  such  as  those  near  1980  cm"1  until  temperature  in  excess 
of  U000°K  are  reached  in  which  case  the  establishment  of  an  unfavorable  chemical 
equilibrium  situation  renders  HgO  essentially  useless  as  a  CO  laser  coupler. 

Deuteroxide  Radical  -  0D 

0 he m  ic  al  _C  h a r act  eris  ti  cs  _ 

Examination  of  Fig.  7  reveals  that  at  temperatures  in  excess  of  5000°K  in 
an  I^-HgO  mixture  the  concentration  of  hydroxide  radical,  OH,  is  comparable  to  that 
of  H20.  The  fundamental  vibrational  frequency  of  this  molecule  is  3735  cm"1,  too 
far  from  2000  cm"1  to  conceivably  be  of  use  as  a  coupler  at  the  latter  frequency. 
However,  the  corresponding  frequency  for  the  deuteroxide  radical  OD,  is  2721  cm'1, 
raising  the  possibility  that  transitions  lying  in  the  far  P  branch  of  the  transition 
may  be  of  use  for  coupling  at  high  temperatures.  As  will  be  seen  in  the  section 
dealing  with  D20,  the  peak  OD  concentration  in  a  mixture  made  initially  of  D20  and 
Hg  is  more  that  an  order  of  magnitude  lower  than  that  shown  for  OH  in  the  correspond¬ 
ing  H20/Hg  mixture  in  Fig.  7  due  to  H-D  exchange  with  the  large  H  atom  population 
present.  For  this  reason  the  most  desirable  chemical  environment  for  the  0D 
molecule  is  a  mixture  made  initially  of  D20  and  D2  and  for  such  a  mixture  the  0D 
concentration  will  be  assumed  to  be  identical  to  the  OH  concentrating  shown  in  Fig.  7 
for  the  corresponding  H20/H2  mixture,  with  the  calculated  loss  per  unit  length  being 
therefore  a  upper  bound. 

Opectros £0£i £  Characteristic^ 

The  ground  state  vibrational  and  rotational  energy  levels  for  0D  were 
calculated  using  the  energy  constants  provided  by  Kerzberg.2  Since  the  ground 
electronic  state  is  a  doublet,  each  substrate  of  which  is  further  split  by  lambda 
doubling,  the  absorption  spectrum  consists  of  groups  of  four  closely  packed  lines  in 
place  of  the  single  lines  which  appeared  in  the  case  of  CO.  In  addition  the  line 
strength  has  a  different  algebraic  form  from  those  of  CO.2^  The  value  of  SBand  for 
OH  is  110  cm-2  amagat"1  3  an<i  this  value  was  scaled  by  the  ratio  of  band  center 
frequencies  of  OH  and  0D  to  81  cm"2  amagat'1.  The  line  broadening  coefficients  were 
taken  from  Ludwig  and  Malkmus^  with  no  further  correction  for  reduced  mass  for 
broadening  by  Dg.  The  value  of  Q(T)  was  calculated  by  summation  over  the  first  60 
rotational  levels  and  first  11  vibrational  levels  of  the  four  electronic  sublevels. 
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convergent  values  of  Q  being  obtained  over  the  entire  range  of  interest,  bine  by 
line  calculation  of  a  for  one  atmosphere  OD  in  10  atmospheres  Dg  was  carried  ut  for 
the  temperature  range  300°K  to  6000°K.  The  calculation  included  P  transitions  up  to 
Pq ( 30 )  for  each  of  the  four  electronic  sublevels,  and  the  results  are  shown  in  Fig. 
14.  At  1000°K  the  three  main  peaks  shown  are  each  made  up  of  transitions  from  each 
the  four  electronic  states  mentioned  aoove  with  the  peaks  centered  on  the  PQ(24), 
P0(25)  and  PQ(2 6)  transitions.  Increase  in  temperature  results  in  significant  con¬ 
tributions  from  excited  vibrational  states.  The  value  of  agp  with  and  without 
equilibrium  chemistry  considerations  included  is  shown  in  Fig.  15  for  the  P3(l3)  CO 
laser  line  at  2012.73  cm-1.  It  is  seen  that  although  OD  is  a  strong  coupler  at  high 
temperatures,  due  in  part  to  the  scarcity,  relative  to  HoO,  of  vibrational  degrees 
of  freedom  which  raise  the  value  of  Q  at  high  temperature,  the  relatively  low  value 
of  agp  at  even  at  its  maximum  value  causes  this  molecule  to  be  of  little  improvement 
over  H2O  as  a  CO  laser  energy  absorber.  It  is  recalled  that  OD  is  to  be  considered 
an  alternate  and  not  a  supplement  to  H2O  as  a  coupling  molecule  due  to  hydrogen  atom 
interchange.  No  significant  contribution  to  the  coupling  to  be  expected  from  D^O 
in  the  D2O  in  the  D2-D2O-OD  system,  for  an  estimation  of  k  indicates  this  frequency 
lies  to  far  in  the  wings  of  the  R  branch  of  the  D20  Vg  transition  to  have  significant 
contributions  to  its  absorption  coefficient  from  the  latter. 

Discussion 


In  this  subsection  three  candidate  molecules  have  been  considered  as  couplers  to 
CO  laser  radiation  in  the  2000  cm-1  region.  Both  line-by-line  and,  where  needed, 
band  model  calculations  have  been  employed  to  estimate  absorption  coefficients,  and 
equilibrium  considerations  have  been  used  to  determine  the  temperature  dependent 
density  of  each.  Of  the  three,  CO  seems  by  far  the  most  desirable  candidate  barring 
the  highly  unlikely  rapid  establishment  of  equilibrium  chemistry  conditions  at 
temperatures  below  1500°K.  This  molecule  exhibits  large  absorption  and  good  chemical 
stability  with  resultingly  good  coupling  (a  >  0.01  cm-1)  properties  up  to  6000°K  and 
in  addition  will  not  exhibit  optical  saturation  in  a  hydrogen  atmosphere  at  the 
temperatures  of  intrest.  In  contrast,  H2O  exhibits  both  a  smaller  absorption  coeffi¬ 
cient  and  poorer  chemical  stability  than  CO,  and  the  D20-D2-0D  system,  conceived 
originally  to  take  advantage  of  the  thermal  decomposition  of  D20,  does  not  equal  CO 
as  a  coupler  because  the  0D  concentration,  even  at  its  largest  is  still  too  small 


to  provide  significant  coupling  compared  to  CO. 
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Wavelength  Region:  10.6  urn 


DoO,  Ho0  and  HDO 


Chemical  Characteristics 


The  equilibrium  chemical  situation  for  mixtures  of  H2  and  H20  has  already 
been  discussed,  and  Fig.  l6  presents  the  corresponding  data  for  a  mixture  made 
up  initially  of  H2  at  ten  atmospheres  and  D20  at  one  atmosphere.  It  is  seen 
that,  were  hydrogen  atom  interchange  readily  accomplished  at  all  temperatures, 
the  DoO  concentration  would  be  rapidly  reduced  at  all  temperatures,  the  resulting 
mixture  being  made  up  of  H20,  HDO  and  P2O.  In  this  system  the  degree  to  which 
chemical  equilibrium  is  established  is  probably  determined  in  large  part  by  the  rate 

of  the  reaction  in  which  a  hydrogen  atom  extracts  a  deuterium  atom  from  DoO  to  form 

PP 

HD  and  0D.  Using  Dixon  Lewis'  rate  constant  for  the  reaction  along  with  the  cal¬ 
culated  hydrogen  atom  concentration  due  to  the  dissociation  as  shown  in  Fig.  l6  it 
is  determined  that  reaction  time  becomes  10  msec,  equal  to  the  residence  time  of 
the  mixture  in  the  high  temperature  region,  at  about  1500°C.  As  a  result  it  is 
assumed  that  a  the  mixture  made  initially  of  1  atm  D20  in  10  atm  H2,  the  D2O  pressure 
will  remain  1  atm  until  1500°C  at  which  temperature  the  calculated  equilibrium 
pressures  will  be  establi shed. 

Li  ne^&yy-  L i  ne^Calculati  ons 

HoO  -  Since  for  all  three  species  of  water  molecule  the  dominant  absorption 
lines  arise  from  either  the  P  branch  of  the  v2  transition  or  from  pure  rotational 
transitions  (type  b),  the  energy  levels  of  interest  in  the  case  of  H20  lie  in  the 
(000),  (010)  and  (020)  vibrational  levels,  and  the  extent  of  the  data  available  for 
this  levels  has  already  been  noted.  For  the  P  branch  of  the  v2  transition,  pre¬ 
cisely  the  same  procedure  was  followed  for  the  frequency  internal  from  895  cm"'*’ 
to  995  cm"*-  as  was  described  for  the  interval  from  1965  cm-^  to  2065  cm-1  so 
that  no  more  need  be  mentioned  here  concerning  the  transitions  except  to  note 
that  Wacker  and  Pratto's^O  calculated  line  strength  are  underestimates  of  the 
true  value  in  the  P  branch  region.  For  pure  rotational  transitions,  only  transi¬ 
tions  within  the  (000)  vibrational  level  were  considered.  For  transitions  not 
included  in  McClatchey's  compendium, ^  the  value  of  MQ,  u>^,  Ej  and  Sjii,  as  cal- 
ulated  from  Wacker  and  Pratto's  tables  and  divided  by  2,  were  substituted  in 

3  to  give  S(296).  The  division  by  2  came  as  a  result  of  trial  calculations  of 
.5(296)  which  were  compared  to  McClatchey's^  results  and  found  to  be  too  high, 
consistent  with  the  expected  effect  of  centrifugal  distortion  on  R  branch  type  line 
strenghts.  As  before,  the  nuclear  spin  statistical  factors  of  1/U  and  3/*i  were 
correctly  applied  in  using  Rq .  (3). 
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In  Fig.  IT,  the  calculated  value  of  a  is  given  as  a  function  of  temperature  in 
the  frequency  range  from  895  cm-^  to  995  cm”'  for  a  mixture  of  one  atmosphere  H^O 
in  ten  atmospheres  H2  neglecting  chemistry  effects.  In  this  frequency  interval,  many 
of  the  most  intense  absorption  lines  arise  from  pure  rotational  transitions  with  the 
result  that  in  the  frequency  region  of  interest  the  peak  intensities  show  a  definite 
decrease  with  increasing  w  similar  to  the  R  branch  behavior  seen  for  H^O  at  2000  cm-1. 
The  decrease  is  not  so  dramatic  as  at  the  larger  frequency  however  due  to  the  presence 
of  v2  r  branch  transitions  in  the  same  frequency  region. 

Laser  transitions  capable  of  furnishing  high  power  levels  and  in  addition  showing 
relatively  low  coupling  to  H20  at  300°K,  and  therefore  good  atmospheric  propagation 
properties,  and  high  coupling  to  H^O  at  higher  temperatures  are  the  P(20)  line  at 
9l*U.15  cm--1  and  the  R(l6)  line  at  972.21*  cm-1-  Good  high  temperature  coupling  and 
good  atmpsheric  propagation  are  indicated  at  962  cm  ,  but  the  corresponding  C02 
laser  line,  the  R(2)  is  not  a  high  power  line. 

HDO  -  As  was  the  case  for  Ho0  the  vibrational  energy  levels  of  interest  are 
the  (000),  (010)  and  (020).  The  (000)  level  is  complete  through  J  equals  11  with 
some  energy  levels  known  for  J  as  high  as  19.^9*  ^ >  7  For  the  (010)  level 
all  levels  for  J  up  to  nine  are  known  with  some  levels  known  for  J  as  high  as  16.^ 
Finally,  for  the  (020)  all  levels  for  J  up  to  four  with  some  levels  known  for 
J  as  high  as  6.^9  Unlike  H20  (or  D20)  HDO  possesses  no  axis  of  symmetry,  a 
fact  which  has  two  results:  first  there  is  no  nuclear  spin  factor  to  be  introduced 
into  the  line  intensity  formula.  Second,  both  the  v2  and  rotational  band 
exhibit  both  type  a  and  type  b  transitions  with  the  former  being  approximately 
half  as  strong  as  the  former.  Thus  the  HDO  absorption  spectrum  exhibits  more 
lines  than  that  of  H20  and  these  lines  are  less  intense.  Due  to  D  being  heavier 
than  H,  all  rotational  and  v2  lines  are  displaced  to  lower  frequencies  compared 
to  their  counterparts  in  h'20. 

The  HDO  line  by  line  calculation  was  carried  out  in  the  same  manner  as  was 
done  for  H20  with  Sg  being  scaled  down  from  its  H20  value  by  the  ratio  of  the 
u>0  values  of  the  v2  band  for  H20  and  HDO.^  The  results  of  the  calculation  are 
shown  in  Fig.  18  as  a  function  of  temperature ,  and  the  properties  described 
above  are  apparent.  The  HDO  spectrum  appears  to  have  more  numerous  but  less 
intense  spectral  lines  than  does  Ho0  and  the  decrease  in  intensity  is  with 
decreasing  frequency,  reflecting  the  greater  dominance  of  v2  transitions  in 
this  frequency  region  for  the  molecule. 

For  the  gas  mixture  of  interest,  namely  D20  at  one  atmosphere  in  10 
atmospheres  H2,  HDO  is  an  intermediate  product  between  the  low  temperature 
coupler,  D20,  and  the  high  temperature  one,  H20.  Thus  no  dovetailing  of  laser 
transition  to  coupler  transition  is  attempted. 
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ABSORPTION  SPECTRA-ONE  ATM  H20  IN  10  ATM  H2 


FIG.  17 
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ABSORPTION  SPECTRA-ONE  ATM  HDO  IN  10  ATM  *2 


FIG.  18 
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^2—  ~  For  the  D2O  molecule,  only  v2  transitions  make  significant  contributions 
to  the  absorption  spectrum  in  the  frequency  range  between  985  and  995  cm-1.  This 
is  the  result  of  the  replacing  of  the  second  hydrogen  atom  with  the  more  heavy 
deuteron  which  decreases  w0  for  v2  to  approximately  ll80  cm--*-  and  shifts  all 
rotational  transitions  to  still  lower  frequencies  compared  to  HgO.  Although  the 
transitions  in  the  manifolds  (OOl)-(Oll),  (011)-(02l)  and  (lOl)-(lll)  were  included 
in  the  D20  calculation, 29  the  dominant  transition  still  involved  the  (OOO)-(OIO)  and 
(010)-(020)  manifolds.  For  the  (000)  level  all  levels  through  J  equals  15  have 
known  values, *-5 »  31  for  the  (010)  level,  all  levels  through  J  equals  6  are  known  and 
some  levels  for  J  as  large  as  12  are  known. 32  No  compilation  of  rotational  levels 
for  (020)  was  found.  In  the  calculations  to  be  presented  first,  the  rotational 
levels  of  the  (020)  level  were  estimated  using  a  method  derived  from  the  observation 
that  in  the  v^*1  vibrational  level  of  a  diatomic  molecule  the  energy  separation  between 
the  Jth  and  zeroeth  rotational  level  decreases  linearly  with  v.  Specifically,  the 
quantity  AEqVq  (J,  t)  was  defined  as  follows 

*Eovo<^  ^  0 V  0  ^  —  ^  0 VO  (0,0)  (18) 

where  E^v0(0,0)  was  calculated  from  the  formula  of  Benedict  Gailar  and  Plyler.^9 
The  energy  level  F02Q(J,  t)  was  then  calculated  according  to 

e020(J,t)  =  E020(0,0)  +  AE0io(J,t)  +(AEoio(J,t)-  AE000(j,r))  (l9) 

In  this  way  rotational  energy  levels  in  the  (020)  level  were  estimated  for  the 
same  (J,  t)  as  are  known  for  the  (010)  level. 

D20  having  an  axis  of  symmetry,  nuclear  spin  factors  again  occur  in  the  expression 
for  the  line  intensity.  For  levels  of  the  type  (++)  or  ( — )  the  factor  is  2/3  while 
for  levels  of  the  type  (-+)  or  (+- )  the  factor  is  1/3. ^  As  in  the  case  of  1I20,  only 
type  b  transitions  occur. 

The  calculation  of  S(296)  was  done  in  the  same  manner  as  for  H20  and  HD0  using 
the  appropriate  DpO  energy  constants  constants  and  scaling  Sg  according  to  w0. 

The  line  by  line  calculation  of  a  using  the  energy  levels  described  above  was 
carried  out  as  a  function  of  temperature,  and  the  results  of  this  calculation  are 
shown  in  Fig.  19.  As  in  this  case  of  HD0,  the  dominance  of  the  \>2  transition,  with 
origin  near  ll80,  is  evident  as  a  exhibits  a  general  decrease  with  decreasing  fre¬ 
quency  typical  of  a  P  branch  wing.  In  addition,  at  the  lowest  temperatures  a  for 
D20  is  seen  to  be  significantly  larger  than  that  for  H2O.  Particularly  interesting 
C02  laser  transitions  from  the  standpoint  of  coupling  are  the  P(l8)  at  9U5  cm-1, 
the  R(28)  at  980  cm-^  and  the  R ( 8 )  at  968  cm”*-. 
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In  contrast  to  HpO,  the  calculated  DpO  spectrum  exhibits  relatively  little 
growth  in  intensity  and  line  density  with  increasing  temperature  in  spite  of  the 
inclusion  of  some  vp  transitions  originating  from  the  (001),  (Oil)  and  (101) 
vibrational  levels.  It  was  suspected  that  this  may  be  due  to  the  relatively  small 
size  of  the  (010)  energy  level  data  base  which  resulted  in  an  inability  to  calculate 
vp  transitions  originating  from  the  (000)  level,  and  a  study  was  carried  out  to  see 
whether  it  was  possible  to  estimate  the  frequency  positions  of  important  high 
temperature  DpO  absorption  transitions  and  thereby  to  include  these  transitions  in 
the  calculation  of  a.  Because  of  the  symmetry  selection  rules  described  previously, 
transition  between  rotational  levels  of  DpO  follow  certain  sequences  which  can  be 
conveniently  denoted  by  (AJ,  At,  J"  +  t")  where  AJ  is  simply  the  difference  J'-J", 
and  At  is  t'-t"  where  the  single  prime  superscript  refers  to  the  higher  lying 
energy  level  and  the  double  prime  to  the  lower  lying  energy  level  involved  in  the 
transition.  Numerical  experimentation  using  the  known  D2O  energy  levels  revealed 
that  the  following  sequences  contained  transitions  in  the  frequency  region  from 
895  cm"1  to  995  cm”1:  (-1,  -6,  8),  (-1,  -6,  10),  (-1,  -U ,  5),  (-1,  -U,  -7),  (-1, 

-1*,  -9),  (-1,  -2,  2n)  where  n  ranged  from  two  to  ten,  (-1,  0,  2n  +l),  where  n  ranged 
from  four  to  nine,  and  (0,  -6,  n)  where  n  ranged  from  6  to  lU.  In  addition  it  was 
found  that  many  of  these  transitions  had  not  been  included  in  the  calculation  of  a 
because  of  lack  of  knowledge  of  the  value  of  Eoio(J’ »  T' The  frequency  of  the 
missing  transitions  was  estimated  by  postulating  that  the  energy  difference  between 

(J',  t')  and  Eq10(J',  -J' )  was  the  same  in  the  (000)  state.  Since  Eq10(J',  -J') 
was  known  for  J  as  large  as  12,  this  method  of  estimation  allowed  the  generation  of 
a  considerable  number  of  transitions.  Employing  these  estimated  energy  levels  and 
Eq.  (19)  permitted  calculation  of  transition  frequencies  between  the  (010)  and 
(020)  levels.  It  must  be  noted  that  the  method  of  estimation  used  is  not  very 
accurate.  For  instance,  the  value  of  Eq1Q(7,2)  calculated  in  this  way  was  1751  cm-1 
compared  to  the  known  value  of  1773  cm-1.  The  point  to  keep  in  mind  is  that  the 
higher  temperature  DpO  spectra  shown  in  Fig.  19  are  missing  many  important  transitions, 
and  the  point  of  this  exercise  is  to  demonstrate  the  effect  of  transitions  from 
higher  energy  levels  in  the  (000)  and  (010)  states  upon  the  absorption  spectrum. 
Finally,  it  is  noted  that  At  values  less  than  -6  were  not  included  due  to  the  small, 
less  than  10_!+,  value  of  their  line  strengths  as  calculated  by  Wacker  and  Pratto.20 

The  D2O  absorption  spectrum,  calculated  with  the  transitions  estimated  as 
described  above  included,  is  shown  in  Fig.  20  as  a  function  of  temperature.  Although 
it  is  again  emphasized  that  these  results  are  not  accurate  detailed  representations, 
the  observed  qualitative  effect  of  expanding  the  D2O  energy  level  data  base  is  most 
likely  real  in  nature,  with  the  spectrum  being  richer  in  lines  as  well  as  being  more 
intense  with  the  intensity  varying  with  temperature  in  a  manner  similar  to  that 
shown  for  H20  earlier.  Reference  to  the  results  in  Fig.  20  will  be  made  in  the 
section  concerning  band  model  calculations  for  D2O. 
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Band  Model  Calculations 


Ho0_  -  in  the  spectral  region  near  1000  cm-l  the  band  model  analysis  must  be 
carried  out  for  two  transitions,  the  v2  band  origin  near  1600  cm-1  and  the 

pure  rotational  band.  Ludwig1*  provides  the  value  of  Ic  for  the  frequency  range  up  to 
1600  cm-1  for  temperatures  as  high  as  3000°K  and  this  data  was  transformed  to  the 
functions  f r (u> )  and  g(iii),  given  by  Eqs .  (9)  and  (10)  in  the  frequency  regions 
below  1000  cm-1  and  above  1000  cm-1  respectively.  For  each  function  the  value  of 
aj*  was  obtained  from  the  position  of  the  maximum  for  each  function  and  the  quantities 
C'  and  C,  where 


87T  3  M02 

3  he 


(20) 


were  obtained  from  the  values  of  g(w*)  and  f(u*)  respectively. 


The  derived  values  at  C  and  w*  for  the  v2  P  branch  are  shown  in  Fig.  21.  Com¬ 
parison  with  Fig.  11  reveals  a  marked  consistancy  between  the  u*  values  for  the  P 
and  R  branches.  In  each  case  the  temperature  dependence,  T1^2,  as  well  as  the 
magnitude  of  ui*  are  nearly  the  same  as  would  be  expected  from  the  definition  of  u* 
given  by  Eq.  (12).  The  value  of  C  for  the  P  branch,  while  agreeing  with  that 
found  for  the  R  branch  at  lower  temperatures  is  seen  to  increase  faster  with 
temperature  to  become  nearly  an  order  of  magnitude  larger  than  the  R  branch  value 
at  300°K.  This  may  be  attributable  to  the  fact  that  the  P  branches  of  transitions 
arising  from  the  (010)  and  higher  excited  states  are  displaced  to  lower  frequencies 
relative  to  transitions  originating  from  the  (000)  level  with  the  result  that  the 
rate  of  increase  of  C  with  temperature  is  accelerated.  In  addition,  due  to  centri¬ 
fugal  distortion  effects,  the  line  intensities  for  transitions  in  the  P  branch 
originating  from  states  having  large  angular  momentum  will  be  higher  than  in  the 
R  branch.  To  the  extent  that  these  transitions  become  increasingly  important  at 
high  temperature,  the  observed  value  of  C  will  be  augmented  over  its  value  in  the 
absence  of  centrifugal  distortion.  Each  curve  is  shown  extrapolated  to  6000°K  and 
the  extrapolated  values  are  used  in  the  band  model  calculation  of  k. 

Figure  22  presents  the  temperature  variation  of  C'  an  for  the  pure  rotational 
transitions  of  H20.  As  expected,  the  value  of  u*  scales  as  T1/2  and  is  essentially 
identical  to  u*  for  the  vp  p  branch  as  expected  from  the  definition  of  w*.  Also 
shown  in  Fig.  22  as  a  dotted  line  is  the  value  of  u*  calculated  from  Eq.  (12) 
using  for  B  the  average  of  the  three  rotational  constants  of  H20,29  17.1  cm-^.  It 
is  seen  that  the  difference  between  the  derived  and  theoretical  values  of  u>*  is 
about  +30  percent  for  the  three  bands  considered.  In  Fig.  22  the  value  of  C' is 
seen  to  be  essentially  temperature  independent  consistent  with  Eq.  (20)  having 
an  average  magnitude  of  about  1*0  cm"-'-  amagat"^.  In  comparison,  the  calculated  value 
of  C' ,  using  Eq.  (20)  and  the  measured  dipole  moment  of  HpO  is  38  cm-1  amagat-^ 
in  substantial  agreement  with  the  results  of  the  band  model  data.  The  extrapolated 
values  of  C'  and  uj*  at  temperatures  above  3000°K  are  used  in  the  band  model  analysis. 
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The  band  model  analysis  for  H20  was  completed  by  calculating  fr(w*)  and  frv(u) 
at  the  desired  temperature  using  the  appropriate  (<u*,  C')  and  (w*,  C)  values, 
adding  the  two  terms  and  multiplying  by  the  stimulated  emission  factor  (l-exp(-hcw/ 
kT)).  The  results  are  given  for  950  cm-1  in  Fig.  23  along  with  Ludvig's^  values  and 
the  value  of  a  (9^5  cm--*-)/u.  The  agreement  between  the  k  and  a/u  is  quite  good 
until  a  temperature  is  reached  where  absorption  by  molecules  in  untabulated  very 
high  rotational  energy  levels  becomes  dominant  in  calculating  a. 


As  was  the  case  at  2000  cm-*,  the  factor  (l  +  ku/ka)^^  is  essentially  unity 
in  the  present  study  and  so  ku  is  plotted  vs  T  in  Fig.  2b  for  the  cases  in  which  u 
either  includes  or  does  not  include  the  equillibrium  chemistry  effects  relevant  to 
an  initial  mixture  of  10  atm  Hp  and  one  atm  HgO.  Taking  a  ku  value  of  10-^  cm-^ 
as  the  point  at  which  H20  is  a  useful  coupler,  it  is  seen  that  H2O  does  not  "turn  on" 
until  about  1300°K  and  is  "turned  off"  by  chemistry  effects  at  i4800°K.  At  2000  cm--*- 
(Fig.  13)  it  is  recalled  that  the  corresponding  temperatures  were  1*80°K  and  3500°K 
and  that  the  maximum  ku  value  was  about  four  times  smaller  than  that  shown  in  Fig.  2b. 


HDO  -  The  band  model  calculation  for  HDO  consisted  first  of  scaling  the  (C,  w*) 
and  (C’,  w*)  sets  obtained  for  the  HgO  molecule  to  correspond  to  the  HDO  molecule. 

That  is  each  w*  value  was  scaled  according  to  B17"  where  B  is  the  average  of  the 
three  rotational  energy  constants  of  the  molecule.  The  value  of  C'  was  taken  to  be 
the  same  for  each  molecule.  Since  Sg  scales  as  u0,  C  is  expected  to  scale  only  as 
qv(T)  and  thus  to  increase  with  temperature  faster  than  does  H20.  However, 
centrifugal  distortion  effects  are  expected  to  be  smaller  in  HDO  than  in  HpO,  there¬ 
by  causing  C  to  increase  not  as  quickly  with  T.  To  take  the  latter  into  account, 
be  value  of  C  used  for  HDO  was  taken  to  be  the  same  as  that  of  H2O. 

The  value  of  k  at  950  cm-*  and  a/u  at  are  plotted  us  temperature  in  Fig.  25. 
t.e  inadequacy  of  the  line-by-line  calculation  data  was  a  temperature  as  low  as  500°K 
for  the  HDO  molecule  is  evident. 

u20  -  As  in  the  case  of  HDO  no  change  was  made  in  the  magnitude  of  C,  and  the 
:.'-r  band  model  parameters  were  scaled  to  the  D2O  molecule  as  directed  by  Eqs. 

md  (20).  The  values  of  k,  a(9^5)/u  for  the  unexpanded  D2O  data  base  and 
u\)b[j)/u  for  the  expanded  base  are  shown  in  Fig.  26  as  functions  of  temperature.  The 
a-: ditioral  PpO  absorption  lines  generated  by  expanding  the  data  base  are  seen  to 
enhance  the  adequacy  of  the  line-by-line  calculation  as  an  estimation  of  the  absorp¬ 
tion  coefficient. 
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FIG.  24 
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Finally,  Fig.  27  presents  rTu  for  the  one  atmosphere  D20/10  atmosphere  H2 
mixture.  The  two  sets  of  curves,  as  labeled,  include  and  ignore  equilibrium 
chemistry  effects.  Since  chemical  equilibrium  is  not  expected  to  become 
established  in  the  practical  device  under  consideration  until  a  temperature  of 
1500°K  is  reached,  the  use  of  DpO  instead  of  H2O  is  seen  to  result  in  significantly, 
more  than  an  order  of  magnitude  at  temperatures  below  1000°K,  higher  coupling  ’until 
temperatures  in  excess  of  1500°K  are  reached.  At  this  point,  with  the  establishment 
of  chemical  equilibrium,  H2O  becomes  the  dominant  coupler  and  remains  so. 

Ammonia 


Ch  eraic  al_Ch  ar  ac_t£r  i_s_t  i  c_s_ 

A  mixture  consisting  initially  of  one  atmosphere  NH3  and  ten  atmospheres 
H0  exhibits  the  equilibrium  composition  shown  in  Fig.  28.  Taking  again  1500°K  to 
be  the  temperature  at  which  chemical  equilibrium  is  reached  within  the  10  msec 
residence  time  of  the  device,  it  is  readily  seen  that  NH^  is  a  potentially  useful 
coupler  only  at  temperatures  below  1500°K  due  to  decomposition  of  ammonia  to  give 
N2  and  H2- 

Spec_tros£0£i£  Characteristics. 

OO 

Taylor 3  provides  till  the  data  needed  to  calculate  S(296)  for  ammonia  for 
all  transitions  of  importance  near  ambient  conditions.  Also  provided  are  the 
values  at  300°K  of  (.075  cm-^  atmosphere)  whose  temperature  scaling  was  taken 
to  be  T1^2.  The  value  of  Yfjj{-  was  taken  from  Tayler33  and  assumed  to  scale  as  T-1. 
The  absorption  spectrum  calculated  using  the  line  by  line  method  is  shown  in  Fig.  29. 
The  spectrum  shows  NH^  to  be  an  extremely  good  coupler  with  a  spectral  richness 
which  provides  a  wide  option  in  the  selection  of  a  laser  lines  to  which  to  couple. 

Figure  30  presents  o  (91*5. 9I+  cm-!)  for  NH3  as  a  function  of  temperature  for  the 
two  cases  in  which  chemical  effects  either  are  or  are  not  included.  Hie  frequency 
is  the  C02  P(l8)  laser  line.  It  is  no  surprise  to  find  that  the  effectiveness  of 
as  a  coupler  is  determined  solely  by  the  question  concerning  the  temperature 
at  which  chemical  equilibrium  is  readily  established.  Thus,  the  molecule  is 
"turned  on"  essentially  at  ambient  temperature  and  "turned  off"  when  a  temperature 
of  1500 °K  is  reached. 

Discussion 


In  this  section,  three  systems,  H2O-H2,  D2O-H2  and  NH^-Hp  were  considered  as 
couplers  to  the  COp  laser.  Using  chemical  stability  and  an  absorption  length  of 
100  cm  (a  =  0.01  cm-1)  as  the  most  desirable  properties,  the  following  conclusion  can 
be  reached  for  these  systems  in  which  the  coupling  molecule  pressure  is  initially  one 
atmosphere  and  the  Hp  pressure  is  ten  atmospheres.  The  DpO-Hp  system  seems  the  most 
attractive  of  the  three,  exhibiting  a  low  "turn  on"  temperature  due  to  the  high 
absorbtivity  of  D2O  and  the  fairly  good,  to  U800°K,  chemical  stability  of  the  HpO 
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FIG.  27 


ABSORPTION  PER  CM  (M950  CM-^u)  FOR  1  ATM  D20  IN  10  ATM  H2  VS  TEMPERATURE 


(U„„  INCLUDES  CHEMISTRY  EFFECTSI 


ATMOSPHERES 


molecule  which  become:',  the  dominant  coupler  as  a  result  of  hydrogen  atom  exchange 
at  temperatures  above  1500°K.  The  latter  quality  is  also  of  course  possessed  by 
the  H0O-H2  system  which  also  shows  however  a  higher  "turn  on"  temperature  of  130'/ 

The  NH^  molecule  is  "turned  on"  at  all  temperatures  considered  but  its  power  chemical 
stability  reculates  it  to  use  only  as  a  low  temperature,  less  than  1500°K,  coupler. 


Conclusions 

Various  candidate  molecules  have  been  considered  analytically  as  couplers  to 
iT  and  C0o  laser  radiation  when  immersed  at  a  pressure  of  one  atmosphere  in  ten 
atmospheres  Ho.  Of  all  the  systems  considered,  the  CO  molecule  as  a  coupler  to  CO 
laser  radiation  is  the  best  candidate.  This  molecule  exhibits  an  absorption  per 
unit  length  of  magnitude  10-2  cm--1-  at  temperatures  as  low  as  300°K  and  a  chemical 
stability  which  permits  this  magnitude  to  be  maintained  to  temperatures  in  excess 
of  f.0O0°K.  In  addition  the  fact  that  Ho  is  a  good  vibrational  deactivator  of  CO 
< •  1 i mi  nates  the  problem  of  optical  bleaching  of  the  molecule.  The  system  DpO-Hp  as 
a  coupler  to  COp  laser  radiation  is  second  only  to  the  CO-Hp-CO  laser  system.  In 
this  case  the  high  absorbtivity  of  D^O  causes  the  coupler  to  attain  an  absorption 
cf  •■fficient  of  10-2  cm  at  about  U00°K  but  dissociation  of  HpO  at  high  temperature 
'■ausei-  the  loss  per  unit  length  to  go  below  10-2  cm--1-  at  about  U800°K.  The  HgO- 
::  ■  system  was  considered  for  both  the  CO  and  COp  laser.  For  the  former  the  HpO-Hp 
system  was  effective  from  l400°K  to  3800°K  while  for  the  latter  this  range  ran  from 
1T>0°L  to  !j800oK  with  a  peak  absorption  coefficient  of  nearly  ten  times  larger  than 
the  former.  The  concentration  of  the  OD  molecule  was  never  high  enough  even  in  the 
yh-Do  system  to  be  an  effective  coupler,  and  thermal  decomposition  prevented  the 
NT.  ,-Hp  system  from  being  useful  at  temperatures  much  above  1500°K. 
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SECTION  III 
EXPERIMENTAL  METHOD 

Int  roduct ion 

The  purpose  of  the  experimental  segment  of  this  program  consisted  of  measuring 
the  temperature  dependent  absorption  per  cm  of  H2/H2O,  H2/D2O  and  H2/NH3  mixtures 
at  10.6  pm  wavelength,  approximately  ten  atmospheres  pressures  and  at  tempera¬ 
tures  as  high  as  6000*K.  To  attain  this  information,  each  mixture  was  placed 
in  a  high  pressure  gas  cell  in  which  a  hot  plasma,  sustained  by  a  focused 
high  power  CW  CO2  laser  beam,  provided  the  region  of  high  temperature.  This 
cell  was  placed  in  one  arm  of  a  Mach  Zehnder  interferometer  with  a  CO2  laser 
light  source,  operating  at  the  wavelength  of  interest,  and  the  interferometer 
was  modified  to  provide  nearly  simultaneous  determination  of  the  power  loss  and 
phase  shift  of  the  light  traversing  the  ceil  due  to  the  presence  of  the  hot  gas 
inside  it.  Due  to  the  axial  symmetry  in  the  properties  about  the  plasma,  these 
data  were  reduced  to  provide  the  radial  dependence  of  both  the  gas  temperature 
and  Che  absorption  per  cm  and  therefore  the  dependence  of  the  latter  on  the 
f  o rme  r . 

As  originally  planned,  measurements  were  also  to  be  taken  at  5pM  in  a  similar 
manner  using  a  CO  laser  instead  of  the  C02  laser  as  the  laser  light  source. 

Although  time  did  not  permit  carrying  out  these  experiments,  the  CO  laser  to  be 
used  was  constructed,  and  this  laser  is  described  in  Appendix  A. 

In  this  section,  we  describe  in  detail  the  experimental  apparatus,  the 
collection  of  data  and  its  reduction  to  provide  the  information  of  interest. 

The  High  Pressure  Cell 

The  gas  cell  used  in  this  study  is  drawn  to  scale  in  Fig.  31.  The  cell  was 
made  from  brass  with  copper  cooling  lines  soldered  to  it.  To  heat  the  cell,  heat¬ 
ing  tape  was  wrapped  around  it.  Between  different  cell  sections  and  at  window-to- 
cell  seals,  0-rings  made  from  a  low  porosity-high  thermal  stability  elastomer 
were  used.  It  is  seen  that  the  cell  is  somewhat  "top  heavy"  in  appearance  due  to 
the  7.5  inch  long  section  which  rests  on  the  cross  which  originally  constituted 
the  cell.  This  section  was  added  to  the  cell  after  an  early  experiment  in  which 
bouyancy  driven  convective  heat  flow  from  the  laser  sustained  plasma  caused  the 
high  power  laser  beam  entrance  window  to  crack.  No  recurrence  of  this  accident 
happened  once  the  entrance  window  was  moved  further  from  the  plasma.  The 
entrance  window  was  zinc  selenide,  two  inches  in  diameter,  0.35  inches  thick  and 
anti-reflectance  ("AR")  coated  for  10.6  pm  wavelength.  The  high  power  beam  was 
directed  through  this  window  by  a  copper  turning  mirror  in  an  adjustable  mount 
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FIG.  31 
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not  shown,  which  was  attached  to  the  top  of  the  cell.  This  beam  had  a  diameter 
of  approximately  1.25  to  1.5  inches  at  this  point,  having  passed  through  a  two- 
to-one  diameter  reducing  telescope  after  emerging  from  the  laser,  the  UTKC 
7.5  kW  high  power  electric  discharge  laser.  The  beam  intensity  was  essentially 
uniform  across  this  diameter  as  the  high  power  laser  oscillates  in  the  so-called 
"top-hat"  mode.  After  passing  through  the  entrance  window,  the  high  power 
laser  beam  is  incident  on  a  one  inch  focal  length  copper  mirror,  located  on 
the  bottom  of  the  cell,  which  retroreflects  the  beam  and  focuses  it  at  a 
point  in  the  center  of  the  field  of  view  of  the  two  side-sets  of  cell  windows. 
Referring  now  to  the  "side  view"  portion  of  Fig.  31,  laser  induced  gas  break¬ 
down  is  made  to  occur  at  the  high  power  laser  focal  point  by  focusing  at  this 
point  the  output  of  a  pulsed  CC^  laser  with  a  four  inch  focal  length  germanium 
lens  through  an  ar  coated  0.35  inch  thick  germanium  window.  With  the  high  power 
CW  laser  beam  quality  optimized  and  this  beam  well  aligned  through  the  input 
window  onto  the  copper  focusing  mirror,  and  with  the  pulsed  breakdown  plasma 
located  near  or  on  the  CW  laser  focal  point,  as  judged  by  visual  observation 
through  the  3/8  inch  thick  pyrex  viewing  window  indicated,  sustaining  of  the 
plasma  with  the  cw  beam  was  quite  readily  accomplished  at  hydrogen  pressures 
from  ten  to  as  low  as  five  atmospheres  so  long  as  the  cell  temperature  was 
near  to  ambient  prior  to  turning  on  the  CW  laser  beam.  At  higher  cell 
temperatures,  obtaining  cw  laser  sustained  plasma  was  much  more  difficult 
regardless  of  the  pressure.  The  reason  for  this  is  not  known  but  may  be  due 
to  the  existence,  at  high  cell  temperatures,  of  temperature  and  therefore 
refractive  index  gradients  in  the  high  power  beam  input  window  with  a  resulting 
decrease  in  the  focusability  of  the  CW  beam. 

The  final  pair  of  windows,  seen  in  the  "front  view"  of  Fig.  31,  are  used  for 
passing  the  probe  laser  beam,  used  for  measuring  the  properties  of  the  gas  mix¬ 
ture,  through  the  cell.  These  windows  are  0.35  inch  thick  germanium,  AR  coated 
for  maximum  transmittance  at  10.6  pm  and  s  polarization  since  the  CO2  laser 
source  is  itself  polarized.  These  windows  are  seen  to  be  displaced  5.6"  from 
the  plasma  with  4"  long  extensions.  The  purpose  of  doing  so  was  to  minimize 
the  degree  of  window  heating  by  the  plasma  which  would  cause  changes  in  their 
absorbi t i vi ty  and/or  refractive  index  during  the  course  of  an  experiment, 
thereby  reducing  the  quality  of  the  data. 


The  Gas  Handling  System 

The  system  used  for  introducing  the  mixture  of  interest  into  the  cell  and 
then  of  disposing  of  it  is  shown  in  Fig.  32.  In  this  figure,  the  gas  flow  pro¬ 
ceeds  from  the  top  of  the  figure  to  the  bottom.  The  source  gases,  argon,  hydro¬ 
gen,  and  ammonia,  are  stored  in  tanks  as  shown.  Argon  was  used  first  to  flush 
air  from  the  cell  and  then  to  raise  the  cell  pressure  to  the  desired  level  prior 
to  the  introduction  of  hydrogen.  The  flow  rate  of  each  gas  was  monitored  using 
the  flow  meter  unit  with  the  pressure  and  species  correction  appropriate  to  each 
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flow  tube  used.  By  monitoring  the  H2  flow  meter  and  the  pressure  gauge  and 
by  adjusting  the  flow  meter  control  valve  the  desired  total  pressure  and  flow 
rate  (10  standard  cubic  feet  per  hour)  were  obtained.  All  measurements  were 
made  under  flowing  conditions. 

For  H2/H2O  (D2O)  experiments  water  was  introduced  to  the  hydrogen  flow 
by  forcing  the  latter  through  the  thermostated  brass  flow  saturater  shown  in  the 
figure.  In  the  saturater,  the  gas  flow  is  forced  through  a  fritted  gas  disper¬ 
sion  tube  immersed  in  the  water  to  form  small  bubbles  which  become  saturated 
with  water  vapor  at  the  vapor  pressure  associated  with  the  saturater  temperature. 

For  H2/NH3  experiments,  ammonia  was  introduced  to  the  stream  by  cracking 
its  flow  meter  control  valve  wide  enough  so  that  the  flow  meter  indicated  that 
NHj  was  flowing  into  the  hydrogen  stream.  For  this  gas,  measurements  were  re¬ 
stricted  to  approximately  5  atmospheres  cell  pressure  since  the  cylinder  delivery 
pressure  was  not  much  higher  than  this  pressure.  Exact  determination  of  the 
NHj  concentration  in  the  flow  was  made  by  flowing  a  known  volume  of  the 
mixture  through  a  known  amount  of  HC1  solution  which  extracted  the  NHj  from 
the  gas  to  form  NH^Cl.  The  amount  of  NH^Cl  formed  and  thus  the  NH-j 
concentration  in  the  flow  was  determined  elect rochemi cally  (Orion  Research 
ammonia  electrode  model  95-10). 

Other  elements  in  the  gas  handling  system  are  as  follows:  a  shut  off  ball 
valve  to  stop  the  flow  of  combustible  H2  in  case  of  failure  of  a  window  or  some 
other  sealing  element  in  the  system,  a  particle  filter  to  remove  particles  larger 
than  60  ^tm  diameter  from  the  flow,  a  condenser  to  remove  water  from  the  gas 
exiting  the  cell,  a  relief  valve  to  prevent  the  occurrence  of  cell  pressures 
higher  than  190  psia,  a  check  valve  to  prevent  the  migration  of  air  into  the 
system  containing  ambient  pressure  hydrogen  in  case  of  a  window  or  seal  failure, 
and  an  electrically  heated  hot  nickel  wire  to  catalyze  the  combustion  of  waste 
hydrogen  emerging  from  the  system. 


Diagnostic  Apparatus 


The  diagnostic  apparatus  to  be  described  was  designed  to  permit  the  nearly 
simultaneous  determination  ot  the  spatial  dependence  of  the  power  loss  and  opti¬ 
cal  phase  shift  experienced  by  10.6  /^m  light  passing  through  the  gas  mixtures 
of  interest  to  this  study.  The  apparatus  is  shown  schematically  in  Fig.  33. 

The  previously  described  cell  was  placed  as  shown  in  one  arm  of  a  Mach  Zehnder 
interferometer  upon  a  sliding  stage  to  permit  the  center  of  the  plasma  formed 
within  the  cell  to  be  moved  horizontally  along  the  y  axis  with  respect  to  the 
interferometer  light  beam.  The  latter  originates  from  a  9W  CO2  electric  dis¬ 
charge  laser  operating  on  the  P ( 1 8 )  transition  at  946  cm-'  frequency.  The 
beam  from  this  laser  passes  through  the  telescope,  made  up  of  mirrors  Mj  (1.0  m 
radius  of  curvature)  and  M2  (l.47m  radius  of  curvature)  to  be  expanded  in  size 
prior  to  entering  the  interferometer  to  a  diameter  of  about  14  mm.  As  the  laser 
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output  was  polarized,  the  interferometer  beam  splitters  Bj  and  B2  were  coated 
for  50  percent  transmission  at  10.6  pim  wavelength  and  S  wave  polarization.  Upon 
emerging  from  the  interferometer  at  beamsplitter  B2  the  C02  laser  beam  passed 
through  lens  Lj ,  whose  function  will  be  described  shortly,  and  onto  the  rotating 
mirror  MR  which  sweeps  the  beam  across  the  gold  doped  germanium  detector  A, 
equipped  with  a  mask  and  a  0.2mro  pinhole,  whose  output  is  thus  the  y  dependent 
intensity  of  the  beam  emerging  at  B2.  To  ensure  that  in  reducing  the  data 
the  intensity  readings  to  be  compared  refer  to  the  same  position  on  the  y  axis, 
the  beam  from  the  helium-neon  (HeNe)  laser  is  also  directed  onto  MR  and  thence 
to  the  pin  diode  detector  P  whose  output  is  superimposed  on  that  from  A,  giving 
a  data  record  consisting  of  a  short,  negative  voltage  spike  corresponding  to  the 
HeNe  beam  followed  by  the  broader  positive  voltage  curve  which  is  proportional 
to  the  CO 2  beam  intensity.  The  mirror  MR  rotates  at  4000  revolutions  per 
minute  so  that  the  time  interval  between  successive  data  records  is  15  msec. 

Under  conditions  where  absorption  of  energy  from  the  high  power  laser  beam, 
by  either  the  plasma  or  simply  the  gas  mixture  alone,  is  high,  density  grad¬ 
ients  within  the  cell  are  large  enough  to  cause  considerable  refraction  of  the 
CO 2  probe  laser  light.  Uncorrected,  this  refraction  causes  large  errors  in 
determining  the  y  dependent  power  loss  and  phase  angle  shift,  and  the  lens  L^ 
is  placed  as  shown  to  reduce  this  problem.  The  germanium  lens  L|  has  a  20  cm 
focal  length,  is  AR  coated  at  10.6  pi m  wavelength  and  is  positioned  to  image  the 
axis  of  the  high  power  laser  within  the  cell  onto  the  pinhole  of  detector  A. 

As  a  result,  all  refraction  effects  occurring  near  the  high  power  beam  axis  are 
removed  as  far  as  their  effect  on  the  measurement  of  the  beam  intensity 

n  / 

profile  at  the  detector  A  is  concerned. 

In  order  to  obtain  both  power  loss  and  phase  shift  information,  two  shutters, 
herein  referred  to  as  shown  on  Fig.  33  were  used.  Shutter  1  is  a  Princeton 
Applied  Research  Model  125A  operating  in  a  manner  so  as  to  chop  the  beam  going 
through  the  cell  at  a  rate  of  2000  per  minute  so  that  this  beam  is  swept 
across  detector  A  on  every  other  revolution  of  the  mirror.  Shutter  1  is  equipped 
with  a  light  emitting  diode  and  detector  which  operate  to  provide  an  output 
voltage  whenever  shutter  1  is  open.  This  voltage  is  fed  into  a  Tektronics  FG504 
function  generator  the  output  of  which  is  used  to  trigger  shutter  2,  a  Laser 
Precision  CTX534  Radiation  Chopper  in  such  a  manner  that  the  latter  is  locked 
on  to  chopping  the  beam  through  the  reference  side  of  the  interferometer  at 
half  the  rate  that  shutter  1  chops  the  beam  passing  through  the  cell.  As  a 
result,  on  four  successive  rotations  of  MR,  detector  A  receives  the  following 
information;  First,  with  both  shutters  open,  the  interferometer  throughput 
I  ^ ( y ) ;  second,  with  shutter  1  closed  and  shutter  2  open,  the  throughput  of 
the  reference  side  of  the  interferometer  only,  I r ( y )  ;  third,  with  shutter  1 
open  and  shutter  2  closed,  the  throughput  of  the  cell  only,  Ic(y);  fourth, 
with  both  shutters  closed  no  laser  intensity  is  seen  by  the  detector,  but  any 
detectable  spontaneous  emission  occurring  in  the  cell  as  a  result  of  the 
presence  of  plasma  is  recorded  by  the  dectector,  In(y). 
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Assuming  that  no  spontaneous  emission  is  detected  so  that  In(y)  is  zero, 
the  absorption  AP(y)  and  phase  angle  change  A <t>( y )  were  calculated  from  intensity 
profiles  measured  with  plasma  present  combined  with  those  obtained  prior  to  the 
establishment  of  plasma,  denoted  by  the  superscript  zero: 


AF-  (y)  -  In  llr  (y)  /  Ic  (y))  +  fn  (1°  ly)/I?  (y)) 


(?}  ) 


A<£ (y) -  cos  '/  I,  (y)-lr  (y)-lc  iy)  \  -  c os'1  /I°  (y)-i°  (y)-Ic°  (y)  \  (2* ) 

\  2(ir  (y)lciy) )l/2  /  v  2  ( (y) )l/2  / 


If  spontaneous  emission  is  detected,  the  quantities  1  ^ ( y )  ,  1 r ( y )  and  I c ( y )  were 
corrected  for  this  simply  by  subtracting  In(y)  from  each.  It  is  thus  seen  that 
the  quantities  Ap(y)  and  A <p ( y )  were  obtained  from  information  collected  in  two 
45  msec  intervals,  respectively  located  after  and  prior  to  the  establishment 
of  the  laser  sustained  plasma. 


Detector  Signal  Conditioning  and  Storage 

As  can  be  seen  from  the  foregoing  discussion,  complete  data  sets  consisting 
of  I  •  ( y ) ,  I  (y),  Ic(y)  and  In(y)  were  collected  at  the  rate  of  about  17  sec  * 
for  a  total  of  nearly  100  sets  for  a  six  second  long  data  collecting  interval, 
and  it  was  necessary  to  find  a  means  of  storing  this  mass  of  data  until  it 
was  convenient  to  analyze  it.  This  was  done  by  storing  the  data  signal  on 
magnetic  tape. 

The  output  of  detectors  A  and  P  were  fed  into  the  positive  and  negative  input 
channels  respectively  of  a  Tektronix  1A7A  differential  amplifier,  and  the  gain  of 
me  latter  was  adjusted  to  obtain  an  amplifier  output  signal  of  0.10  to  0.15  volt 
at  the  maximum  value  of  I-(y).  Larger  output  signals  from  this  amplifier  re¬ 
sulted  in  a  saturated  signal  being  recorded  on  the  tape  due  to  saturation  of  the 
amplifier  of  the  tape  recorder  which  had  a  gain  of  about  ten.  The  output  of  the 
1  A/A  was  fed  into  a  Tektronix  Type  0  operational  amplifier  with  unity  gain  which 
provided  the  impedance  matching  necessary  to  obtain  successful  recording  of  the 
signal  on  the  magnetic  tape  by  the  Sony  AV-3650  videocorder  used.  If  for  the 
interferometer,  the  angle  ♦  is  defined  according  to  Eq.  23, 

I,  *  I  r  +Ic+2  Mrlc^  COS  <fi 


then  the  largest  signal  recorded  is  I^(y)  corresponding  to  ♦  equaling  zero. 
The  magnitude  of  this  signal  was  on  the  order  of  one  volt. 
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Data  Processing 


Data  Selection 


Experimentally  it  was  found  that  subsequent  to  a  significant  change  in 
the  rate  of  energy  depositon  into  the  cell  contents,  either  by  simply  turning 
on  the  high  power  laser  or  by  establishing  a  CW  laser  sustained  plasma  in  the 
cell,  there  was  always  a  period  of  no  more  than  5  seconds  during  which  the  quan¬ 
tity  I-(y)  exhibited  marked  temporal  variation  prior  to  settling  into  a  quasi 
steady  state  in  which  its  temporal  variation  was  much  smaller.  This  induction 
time  is  consistant  in  magnitude  with  the  calculated  time,  t,  needed  to  obtain 
a  steady  state  temperature  distribution  by  thermal  conduction  inside  a  radius 
of  magnitude  a  in  a  gas  of  density,  P,  specific  heat  Cp,  and  thermal  conduc¬ 
tivity  K  as  given  by 


r  -oC£ LQ 
4k 


2 


(2M 


For  hydrogen  at  eleven  atmospheres  and  a  radius,  1.6  cm,  consistant  with  the 
high  power  laser  beam-to-cell  wall  distance,  a  value  of  5  seconds  is  calculated 
for  t.  Typically,  then,  an  experimental  run  that  lasted  15  seconds  provided  10 
seconds  of  data  during  which  a  quasi  steady  state  situation  persisted  inside  the 
cell.  During  this  period,  I^(y)  did  not  remain  constant  but  typically  showed 
periods  180-235  msec  long  of  rapid  variation  followed  by  periods  of  similar  dur¬ 
ation  during  which  it  remained  nearly  constant.  The  reason  for  this  behavior  is 
not  known  but  may  be  due  to  instabilities  in  the  buoyancy  driven  vertical  gas 
flow  field  about  the  high  power  laser  axis.  The  tape  record  for  the  quasi  steady 
state  time  interval  was  scanned,  advancing  the  tape  manually  with  the  recorder 
at  stop  action  playback,  and  data  was  selected  for  reduction  under  the  criterion 
that  the  data  selected  occured  between  successive  identical  l-(y)  profiles. 

Transfer  of  Data  to  Computer 

Each  data  record,  consisting  of  the  He-Ne  and  COj  probe  laser  intensity 
profiles  was  200usec  long  for  the  beam  dimension,  mirror  MR  rotation  rate  and 
MR-to-A  distances  used  in  these  experiments.  From  each  data  record,  detailed 
and  accurate  intensity  information  had  to  be  obtained  for  use  in  Eqs.  (21)  and 
(22).  To  do  this,  each  selected  200usec  long  voltage  trace  stored  on  the  mag¬ 
netic  tape  was  digitized  and  stored  as  a  512  element  array  on  an  IBM  Diskette 
1  type  "floppy  disk."  With  the  tape  recorder  in  the  stop  action-playback  mode 
a  portion  of  the  tape  containing  the  200usec  interval  of  interest  was  repeat¬ 
edly  scanned  and  fed  into  a  Tektronix  7A22  differential  amplifier  in  a  Tektronix 
7704A  oscilloscope  system  equipped  with  a  P7001  processor.  The  desired  200usec 
interval  was  selected  using  the  oscilloscope's  7B71  Delaying  Time  Base,  and  this 
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interval  was  sampled  twenty  times  to  average  out  noise  not  present  on  the  tape 
record  and  then  stored  on  the  "floppy  disk."  The  computer  system  used  is  a 
Digital  RX01  equipped  to  interface  with  the  P7001  processor,  and  the  computer 
program  ST0R2  used  is  listed  in  Appendix  B.  Details  of  the  language  and  com¬ 
puter  subroutines  used  in  this  and  the  other  computer  codes  listed  in  the 
appendices  are  found  in  Ref.  35. 


Final  Adjustments  on  Computer-Stored  Data 

Prior  to  their  use  in  Eqs.  (21)  and  (22),  the  data  records,  now  consisting 
of  512  element  numerical  arrays,  denoted  herein  simply  by  the  alphabetic  character 
associated  with  the  physical  quantity,  are  finally  adjusted  to  ensure  first  that 
the  given  address  in  each  array  refers  to  the  same  value  of  y,  second  that  the 
base  line  of  each  array  is  correctly  zeroed  and  third  that  a  tape  recorder 
imposed  noise  signal  approximately  7  mv  in  magnitude  and  66hsec  in  period  is 
subtracted  from  each  array.  When  appropriate,  the  Spontaneous  emission  from 
the  cell  as  stored  in  In  was  subtracted  from  the  intensity  record.  This 
exercise  is  carried  out  using  the  computer  program  RED5 ,  described  in  Appendix  C. 

B.  Calculation  of  Absorpt ion 

With  the  data  corrected  for  zero  level,  periodic  noise  and  array  address  and 
spontaneous  emission  as  described  above,  a  power  loss  array  was  calculated 
using  a  computer  program  corresponding  to  Eq.  (21),  and  this  program,  DALI,  is 
described  in  Appendix  D. 


Calculation  of  Phase  Angle  Shift 


To  calculate  the  phase  angle  shift  array  experienced  by  the  probe  laser 
beam  passing  through  the  cell  with  plasma  present,  three  distinct  calculations 
were  made,  starting  first  with  the  calculation  of  the  phase  angle  array  ♦  from 
Eq.  (23),  and  the  computer  program,  P,  associated  with  this  calculation,  is  given 
in  Appendix  E.  The  quantity  ♦  was  then  calculated  using  the  computer  program, 
DRF.DA ,  listed  in  Appendix  F.  Operationally  this  program  first  computes  the 
arrays  X  and  cosX  and  then  compares  each  element  in  cos4>  with  each  in  cosX,  to 
select  the  element  in  the  latter  closest  in  value  to  the  former  and  finally  to 
assign  the  corresponding  element  in  X  to  ♦. 


As  indicated  by  Eq.  22,  to  calculate  the  preceding  operations  were  carried 
out  on  data  sets  taken  prior  to  and  during  the  presence  of  the  laser  supported 
plasma,  with  £♦  being  the  resultant  difference  between  the  calculated  values. 
However,  the  only  values  of  ♦  directly  calculable  lie  between  -n  and  0  and  there 
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is  therefore  a  second  array,  corresponding  to  the  same  values  of  cost. 

There  is  therefore  the  ambiguity  as  to  which  of  the  four  quantities  ,  ♦♦♦*, 

-  $  -  ♦*,  or  -  ♦  ♦  ♦*  accurately  describes  A$. 

The  above  ambiguity  is  resolved  by  considering  the  physical  situation  giving 
rise  to  A$.  The  phase  shift  is  given  in  terms  of  the  change  in  the  index  of 
refraction,  n,  along  the  path  (-XQ,  XQ)  traversed  by  the  beam 


A<*>  ly)  -  2JI  |X° 
^  -Xq 


In  U,y)-n°U,y)  )  dx 


(2?) 


where  \  is  the  radiation  wavelength,  10.6pm.  The  index  of  refraction  is  related 
to  the  gas  pressure,  p.  and  temperature  T  as  well  as  the  polarizability  of  the 
gas  a  according  to 


n  -  2  n  pq 
kT 


(261 


It  is  seen  that  A$  (y)  is  negative  and  if  as  in  the  present  experiments,  the 
temperature  distribution  has  cylindrical  symmetry  with  a  maximum  on  the  axis, 
which  is  the  high  power  laser  beam  axis,  then  A<Ky)  passes  through  a  minimum 
at  y  equals  zero.  With  this  knowledge,  A$  is  obtained  first  by  examining  the 
four  quantities  described  above  and  deciding  which  is  most  likely  to  describe 
A«  based  on  its  shape  and  the  position  of  its  minimum  compared  to  the  position 
of  the  plasma  center  as  determined  from  the  position  of  maximum  detected  plasma 
spontaneous  radiation.  This  procedure  was  carried  out  in  the  computer  program 
CHOOS2  described  in  Appendix  G. 

Obtaining  Absorption  Coefficient  and  Temperature 

Given  a  sample  with  axial  symmetry  in  the  property  £  about  the  z  axis  on  any 
xy  plane,  the  observed  quantity  E  measured  along  the  X  direction  at  a  distance  y 
from  the  zx  plane  is  given  by 


iy>  =  2  f 


(r02-y2 


1/2 


(r)  dx 


(27) 
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where  rQ  is  the  radial  boundary  past  which  the  property  t(r)  vanishes.  With 
E(y)  measured  and  known,  t( r)  is  obtained  from  the  Abel  transform"^ 


«(r)  -  -J_  f'o  dE  (y)  ( y2  -  r  2) 1/2  d y  (28) 

n  -'r  dy 

When  the  quantity  E(y)  is  the  absorption  Ap,  «(r)  is  the  difference  between 
the  product  of  the  absorption  coefficient  k(r)  and  the  absorber  density 
u(r)  from  its  value  prior  to  establishing  the  plasma.  Similarly,  when  E(y) 
is  the  phase  angle  shift,  «(r)  is  the  refractive  index  differential  n(r)-n°(r) 

In  the  present  work,  this  technique,  commonly  referred  to  as  Abel  inversion, 
was  accomplished  by  first  fitting  E(y)  to  the  form 

2n 

OnY  (29) 


N 

E(y)=  I 


n  =  o 


for  which  case 


N 

«(r)=  I  2nanFn  (r) 

Or  | 


(30) 


where 


Fn^r.ro) 


(2n-i) 


r?) 


172 


dy' 


and 


,  .  1  #  .m  2m  ,  .  ,  2  2  21n-m-l) 

Fn(y,r)  =  I  H)  2  m1  (n-i)1  (y  ~r  )  2  y _ 


(31) 


That  Eq.  (31)  is  true  is  seen  by  observing  that 


dFn  (y.r0)  v2n‘'  (32) 
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The  task  of  fitting  AP(y)  and  A>(y)  to  y  is  accomplished  by  use  of  the  computer 
program  MASAG5  which  is  listed  in  Appendix  H.  In  this  program,  the  portion 
of  the  array  to  be  fitted  to  a  polynomial  is  selected,  and  random  noise  in  this 
portion  is  averaged  over.  The  subroutine  assigns  values  of  y  to  the  locations 
of  the  array  to  te  fitted,  and  the  z  axis  is  in  practice  always  assigned 
to  whatever  address  is  observed  to  occur  at  the  plasma  center  as  observed  from 
the  spontaneous  emission  from  a  plasma  established  in  H2 .  In  assigning  y 
values  to  array  locations,  use  is  made  of  a  factor  which  is  the  number  of 
centimeters  per  unit  change  in  array  address  number  and  is  based  on  the  rota¬ 
tion  rate  of  the  mirror,  the  mirror-to-detector  distance  and  the  known  200usec 
temporal  width  of  the  arrays  used  in  this  work.  The  size  of  this  factor  is 
.00535  cm  per  array  address  change.  It  is  thus  seen  that  the  200W  pinhole 
causes  the  recorded  data  to  be  effectively  an  average  over  nearly  four  array 
addresses.  In  selecting  the  expression  E(y)  to  fit  the  data,  the  polynominal 
order,  is  increased  until  the  resulting  graph  of  E(y)  is  judged  by  visual 
inspection  to  be  a  good  representation  of  the  data. 

The  value  of  rQ  is  the  value  of  y  for  which  the  polynomial  fit  first 
crosses  zero.  In  the  case  of  A$,  the  procedure  as  described  thus  far  does  not 
permit  accurate  determination  of  the  magnitude  of  the  minimum  in  A$.  That  is, 
the  minimum  calculated  may  differ  from  the  actual  value  by  an  increment  equal 
to  an  integral  factor  of  211.  For  experiments  in  which  the  high  power  laser 
beam  was  simply  focused  into  the  gas  mixture  with  relatively  small  interaction 
with  it,  the  temporal  variation  in  I^  was  slow  enough  to  permit  accurate 
estimation  of  A$,  and  it  was  found  that  in  these  cases  the  magnitude  of  A$  was 
such  that  rQ  was  equal  to  or  slightly  less  than  the  distance  of  the  z  axis  to 
the  cell  wall,  1.6  cm.  When  the  high  power  beam-mixture  interaction  was 
intense  as  in  the  pressure  of  plasma,  the  variation  in  1^  was  too  fast  to 
follow  from  the  recorded  data  and  so  the  increments  of  -2 n  were  added  to  the 
value  of  the  minimum  in  A$  until  the  calculated  value  of  rQ  fell  close  to  1.6 
cm . 


With  the  data  fitted  to  the  polynomials,  the  task  remained  to  carry  out 
the  Abel  inversion,  and  this  was  done  for  both  Ap(y)  and  A*(y)  by  the  computer 
program  ABELS  listed  in  Appendix  I. 

Obtaining  k(r)u(r)  from  Ap(y)  was  carried  out  in  a  quite  straightforward 
manner,  but  the  method  of  obtaining  T(r)  from  &Wy)  merits  more  discussion  at 
this  point.  Combining  Eqs.  (25)  and  (26)  and  substituting  the  density  p  for 
p/kT  inversion  results  in  calculation  of  the  temperature  itself  but  rather  the 
product  of  the  density  and  the  species  averaged  polarizability  ^  of  the 
medium  at  radius  r: 


(  33) 
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where  c^(r)  is  the  value  of  £(r)  with  E(y)  being  4>  (y).  Using  the  temperature 
dependent  chemical  equilibrium  calculation  shown  in  Fig.  7  and  a  similar  one 
for  ten  atmospheres  pure  along  with  literature  values  of  the  zero  fre¬ 
quency  polarizability  of  I^,  H  and  0^  and  I^O^,  the  species 

averaged  polarizability  was  calculated  for  ten  atmospheres  H2  and  compared  to 
that  for  ten  atmospheres  H9  and  0.5  l^O.  The  latter  was  found  to  exceed 
the  former  by  8  percent  due  to  the  large  h^O  polarizability,  but  this  differ¬ 
ence  vanished  at  around  3000*K  due  to  the  dissociation  of  l^O.  Due  to  the 
small  size  of  this  difference,  the  polarizability  of  a  mixture  originally  pure 
H2  was  used  in  the  temperature  calculation.  The  density  P°  in  Eq.  (37)  was 
obtained  from  readings  of  cell  pressre,  PQ,  and  temperature,  T  ,  using  the 
pressure  gage  in  Fig.  32  and  a  chrome  1 -a  1 nme  1  thermocouple  in  the  lower  thermo¬ 
couple  well  in  Fig.  31.  In  practice,  the  cel)  pressure,  P',  during  the  experi¬ 
mental  run  was  found  to  differ  from  P0  by  as  much  as  15  percent  due  to  the 
heat  addition  from  the  high  power  laser  beam.  Accordingly,  this  pressure  was 
used  in  obtaining  the  temperature  at  r  from  p(r)a(r)  and  a  calculated  plot 
of  P(T)5(T)  vs.  T. 

As  a  result  of  this  inversion,  both  the  temperature  and  the  absorption 
per  cm  were  known  as  a  function  of  radial  position  relative  to  the  high  power 
laser  beam  axis,  thereby  revealing  the  dependence  of  the  latter  upon  the 
former. 
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SECTION  IV 

EXPERIMENTAL  RESULTS 


Refraction  Effects 

In  an  experiment  in  whch  a  laser  beam  is  used  to  obtain  spatially  detailed 
information  near  a  hot  dense  plasma  which  gives  rise  to  steep  temperature  and 
refractive  index  gradients,  the  question  arises  as  to  what  extent  the  latter 
cause  significant  refraction  of  the  diagnostic  laser  beam  and  thereby  affect 
the  quality  of  the  data.  As  mentioned  earlier,  the  20  cm  focal  length  lens  was 
placed  at  the  exit  of  the  interferometer  to  correct  such  refraction  effects, 
and  this  section  describes  the  success  obtained. 

With  the  cell  filled  with  H2  at  150  psig,  a  laser  sustained  plasma  was 
established  and  data  taken  as  described  in  the  previous  sections.  Figure  34 
shows  the  arrays  Ir>  Ic,  Ir“  and  I£*  the  former  pair  exhibiting  substantial 
distortion  due  to  spontaneous  emission  from  the  plasma,  and  Fig.  35  shows  the 
array  I  containing  the  detected  spontaneous  emission.  The  arrays  If  and  lc 
were  processed  through  RED5  which  subtracted  out  In  to  eliminate  the  effect 
of  the  emission  on  the  latter,  and  the  same  four  arrays  are  re-plotted  in  Fig. 
3b .  The  attenuation  of  lr  by  the  plasma  is  seen  to  occur  at  nearly  the  same 
location  as  the  onset  of  the  radiation  intensity.  The  exact  separation 
between  the  two  onsets  is  about  0.05  cm,  with  essentially  no  distortion  in  Ir 
occurring  at  points  further  away  from  this  point  which  is  0.4  cm  from  the 
plasma  core.  When  this  data  was  fully  reduced,  the  onset  of  absorption  was 
calculated  to  occur  at  a  radius  of  about  0.4  cm,  where  the  temperature  was 
calculated  to  be  8200'K,  in  substantial  agreement  with  the  results  given  for 
the  onset  of  inverse  bremsstrahlung  absorption  in  hydrogen  in  Ref.  41.  Of  more 
importance  however,  the  fact  that  no  distortion  of  Ir  occurs  at  greater 
radial  positions  and  lower  temperatures,  as  expected  for  pure  H2 ,  infers  that 
refraction,  even  in  the  presence  of  plasma,  has  no  significant  effect  on  the 
data  in  the  temperature  region  of  interest  to  this  study. 


The  Laser  Sustained  Flame 


Int  roduct ion 


In  all  three  of  the  gas  mixtures  reported  on  here,  an  interesting  pheno¬ 
menon  was  observed  which  has  important  implications  for  the  laser  energized 
rocket  thruster  as  well  as  providing  a  second  source  of  experimental  data  for 
the  present  study.  In  all  of  the  l^/NHj  experimental  runs  and  in  many  of  the 
and  l^/DjO  runs,  a  tenuous  visible  glow  appeared  at  the  focal  spot  of 


74 


DETECTOR  SIGNAL  -  mV 


R 79-922895-25 


7 

/ 


MEASURED  SPONTANEOUS  EM1SSK 


-250 


R79-922895-25 


the  high  power  laser  beam  in  the  absence  of  plasma.  This  glow  was  not  observed 
in  pure  1^  and  its  color  was  characteristic  of  the  mixture  in  each  case, 
having  in  each  case  a  shape  reminiscent  of  the  flame  from  a  match  or  stove 
burner,  hence  the  application  of  the  term  "flame"  to  the  phenomenon.  Each 
flame  was  blunt  at  its  l  se  which  was  several  millimeters  in  diameter  and 
located  at  the  high  power  beam  focal  spot.  In  length,  each  flame  extended 
several  cm  above  its  base.  The  term  flame  does  not  infer  the  presence  of 
combustion  in  the  cell  but  rather  the  presence  of  a  very  hot  region,  as  high  as 
1000* K-4000* K ,  heated  by  the  absorption  of  energy  from  the  focused  high  power 
laser  beam  by  the  mixture  under  study. 

Importance 


The  implication  of  this  phenomenon  to  the  laser  energized  rocket  is  seen 
by  examining  the  absorption  per  cm/temperature  curve  for  t^/l^O  in  Kig.  24. 

It  is  seen  that  the  magnitude  of  the  absorption  per  cm  remains  relatively  low 
until  a  temperature  of  about  1500°K  is  reached,  implying  that  some  method  of 
pre-heating  of  the  mixture  is  necessary  to  render  it  suitably  opaque  for 
practical  use.  The  flame  phenomenon  however,  implies  that  a  300*K  l^/l^O 
mixture  irradiated  with  high  intensity  CC^  laser  radiation  is  unstable,  with 
the  absorption  of  a  small  amount  of  laser  power  causing  elevation  of  tempera¬ 
ture  which  causes  an  increase  in  ku  which  causes  more  power  absorption  and  so 
on  until  temperatures  of  practical  interest  are  reached  without  the  use  of 
auxiliary  preheaters  or  low  temperature  absorbers  such  as  NH-j,  which  in  fact 
was  originally  included  in  this  study  because  of  its  possible  use  in  this 
way . 

Water  Flames 


As  mentioned  previously,  plasma  was  most  easily  sustained  when  the  cell 
was  relatively  cool  and  it  was  under  these  conditions  that  l^/t^O  flames 
were  most  readily  observed.  This  fact  may  be  attributable  at  least  in  part  to 
the  greater  focusability  of  the  high  power  laser  beam  when  the  cell  is  cool, 
as  hypothesized  previously.  Another  explanation  may  be  the  formation  of  small 
droplets  of  water  in  the  mixture  when,  as  was  often  the  case  when  water  partial 
pressures  of  about  0.3  atmospheres  were  desired,  the  saturator  temperature 
.  needed  the  cell  temperature.  Evaporation  of  these  droplets  or  larger  ones  on 
f'lusing  mirror  surface  would  cause  the  1^0  concentration  and  thus  the 
•prion  per  cm  in  the  beam  to  increase  substantially  and  create  a  flame. 

•  •  mperature  elevation  characteristic  of  a  flame  was  obtained  in  an 

•f  uhi  1  the  cell  temperature  exceeded  that  of  the  saturator, 

,  .•  «  •  •  m»  o.l  1 1  i  -ns  such  as  the  presence  of  droplets  may  not  be 

•  fiat  1  n  Experiments  with  and  Ar/l^O  mixtures 

'  w*.  is  )er. -rnuned  bv  the  majority  host  gas,  for 
*  . .  '  t  .  «r  r  -n  h  t  '  s,  a  white  center 
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The  l^/NH-j,  flame,  like  that  for  H2/H2O  and  H2/D20  was  orange  in  color  but, 
for  the  NH^  concentration  tested,  redder  and  not  as  bright. 


Example  of  Data  Reduction 


Reduction  of  Data 


With  the  cell  filled  with  hydrogen  to  a  pressure  Pp  of  147  psig  (11.0  atmos¬ 
pheres)  H2  along  with  0.3  atm  1^0,  the  recorded  detector  output  appeared  as  shown 
in  Fig.  37,  and  the  presence  of  plasma  in  the  cell  had  the  effect  shown  in  Fig. 
38.  Although  the  detected  spontaneous  emission  in  this  experiment  was  not 
large  enough  to  permit  determination  of  the  position  of  the  high  power  laser 
beam  axis,  this  property  can  be  obtained  from  the  calculated  values  of  A4>  and 
AP  shown  in  Fig.  39  and  is  seen  to  be  at  array  address  170. 

Part  of  the  process  of  fitting  an  analytic  expression  such  as  that  given 
in  Eq.  33  to  the  A<|>  or  AP  data  is  choosing  the  array  address  interval  over 
which  the  fitting  is  to  be  made.  The  main  criterion  for  this  selection  process 
is  that  the  detector  s i gna 1 -t o-noi se  be  large  enough  that  the  resulting  random 
noise  in  A$  and  AP  be  small.  The  resulting  fits  are  shown  in  Figs.  40  and  41, 
and  the  corresponding  curve  of  ku  vs  T,  which  resulted  from  inverting  these 
fits  to  get  the  radial  dependence  of  ku  and  T,  is  marked  with  an  asterisk  in 
Fig.  42 . 

Error  Anal y  s  i  s 


The  inversion  of  the  above  data  was  repeated  with  small  variations  in  each 
of  the  following  measured  input  variables:  AP,  A$ ,  Tq  and  PQ  and  P'.  The 
resulting  change  in  the  calculated  values  of  ku  or  T  reflected  the  sensitivity 
of  these  on  the  variables. 

The  effect  of  changing  AP  is  to  change  the  value  of  rQ  for  AP.  The  applied 
change  in  AP  was  +0.1  which  changed  rQ  by  +.15  mm.  The  resulting  change 
in  ku  was  +0.05  cm-1  on  the  high  power  beam  axis  where  the  temperature 
was  highest.  Far  from  the  axis,  the  variation  in  ku  increased  markedly,  being 
about  +0.15  cm-1  at  the  largest  radii  for  which  the  caluclated  ku  value  exceed 
zero.  Therefore,  the  effect  on  the  ku  vs  T  curve  in  Fig.  42  of  increasing  AP  by  0.1 
is  to  decrease  significantly  the  severity  of  downward  trend  in  ku  with  decreasing 
temperature  at  the  lowest  temperatures  while  having  much  less  effect  at  higher 
t empe ra t ures .  Similarly  decreasing  AP ,  increases  this  downward  tread. 
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Similarly,  the  effect  of  shifting  is  to  change  ru  for  A$  At 
large  radii  and  the  low  end  of  the  temperature  range,  near  2000*K,  it  was  found 
that  a  ±1%  change  in  A$  resulted  in  a  +  3%  change  in  T  while  at  the 
smaller  radii  and  larger  temperatures  the  error  increased  to  +72.  From 
Fig.  41,  it  is  seen  that  the  error  in  At  is  likely  to  be  typically  +3% 
generating  a  temperature  error  of  180°K  near  2000*K  and  about  1200*11  near 
6000°K . 

A  variation  of  +  1K*  in  the  value  of  TQ  as  measured  by  a  thermocouple 
in  the  lower  thermocouple  well  in  Fig.  31  was  calculated  to  generate  a  +1% 
error  in  T  at  2000°K  and  a  ±5%  error  at  6000*K.  Due  to  the  large  thermal 
conductivity  of  hydrogen,  the  massiveness  of  the  cell  and  the  depth  of  the 
well,  it  is  felt  that  the  error  in  TQ  seldom  exceeded  +  1°K. 

A  variation  of  +  1%  in  the  pressure  PQ  and  P'  as  measured  by  the 
pressure  gauge  shown  in  Fig.  32  generated  a  +4%  error  in  T  near  2000*K  and 
a  +12%  error  near  6000°K,  with  the  bulk  of  this  error  being  due  to  the  PQ 
error  rather  than  P1.  The  random  variation  of  the  latter  during  an  experimental 
run  was  observed  to  be  +2  psi  at  11  atmospheres  or  +1.2%  which  in  turn 
generated  +1.2%  in  temperature  over  the  entire  temperature  range.  Of  more 
significance  is  the  systematic  error  in  the  PQ  reading  which  for  this  type 
of  gauge  is  typically  as  high  as  +3%.  The  largest  systematic  error  in  T  is 
in  this  case  +12%  (240*K)  at  2000*K.  and  +36%  (2200*K)  at  6000°K  with 
this  error  being  systematic  rather  than  random  in  nature.  Qualitatively,  using 
the  ku  vs  T  curve  shown  in  Fig.  24  for  comparison,  if  the  pressure  gauge  error 
is  positive,  the  resulting  ku(T)  curve  will  have  a  maximum  at  an  erroneously 
low  temperature  and  too  small  a  half  width  while  a  negative  pressure  error 
will  shift  the  maximum  to  higher  T  and  increase  the  half  width. 


Measured  ku(T)  Curves 


Results  for  H20 


As  mentioned  previously,  in  all  the  experiments  involving  water,  the 
latter  was  added  to  hydrogen  by  bubbling  it  through  a  thermostated  reservoir 
filled  with  the  liquid  prior  to  being  flowed  into  the  cell.  The  cell  was 
generally  kept  near  room  temperature  as  plasma  sustainment  was  found  to  be  most 
attainable  under  this  condition.  However,  under  this  condition  the  condensa¬ 
tion  of  water  inside  the  cell  became  a  problem  which  was  circumvented  by 
imposing  a  limit  on  the  time  for  which  water  was  added  to  the  H2  prior  to 
establishing  plasma  and  taking  data.  This  time  limit,  along  with  the  condensa¬ 
tion  would  act  to  cause  the  indicated  water  pressure,  which  corresponds  to  the 
vapor  pressure  at  the  reservoir  temperature,  to  be  an  upper  bound  on  the 
pressure  actually  in  the  cell.  On  the  other  hand,  to  the  extent  that  H20. 
condenses  on  the  focusing  mirror  and  is  flash  evaporated  into  the  gas  mixture 
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when  the  high  power  iaser  is  turned  on,  the  indicated  H.,0  pressure  would  be 
lower  than  that  actually  in  the  cell.  Thus,  the  actual  water  pressure  in  the 
experiments  reported  on  in  this  and  the  following  D2O  section  is  characterized 
as  being  only  an  estimate. 

The  measured  value  of  ku  as  a  function  of  temperature  is  presented  for  the 
system  in  Fig.  42.  In  this  figure,  there  are  five  separate  curves. 

The  three  solid  curves  are  experimental  results  from  data  obtained  in  the 
presence  of  the  laser  sustained  plasma,  the  one  marked  with  an  asterisk  having 
been  obtained  as  described  in  the  preceding  subsection.  The  dashed  curve  is 
also  an  experimental  result  obtained  in  the  presence  of  the  laser  sustained 
H2/H2O  flame.  The  dotted  curve  is  the  analytical  result  shown  in  Fig.  24. 

Comparing  the  two  sets  of  experimental  results,  denoted  "plasma  data"  and 
"flame  data",  it  is  seen  that  they  are  quite  consistant  in  magnitude  while 
being  complementary  with  respect  to  temperature,  combining  as  they  do  to  cover 
the  temperature  range  from  700°K  to  5500°K. 

Optical  Saturation  Effects 

The  similarity  in  magnitude  between  the  flame  and  plasma  data  at  3000°K  is 
interesting  when  one  notes  that  in  the  former  case,  this  temperature  is  highest 
in  the  cell  and  occurs  at  the  highest  incident  high  power  laser  intensity,  one 

O 

to  ten  MW/cm  ,  whereas  this  same  temperature  range  m  the  case  of  the  plasma 

data  lies  outside  the  core  of  the  laser  sustained  plasma  where  the  incident 

9 

high  power  laser  intensity  is  on  the  order  of  0.01  MW/cm  .  It  is  recalled 
that  the  high  power  laser  operates  at  essentially  the  same  wavelength  as  the 

diagnostic  laser  so  that  the  measured  value  of  ku  is  apparently  not 

greatly  affected  by  optical  saturation  effects  at  intensities  as  high  as  1  to 

10  MW/cm2. 

With  respect  to  optical  saturation,  two  subject  areas  have  been  considered 
during  the  course  of  this  study:  first,  collisional  effects  wherein  the  depopu¬ 
lation  of  an  absorbing  energy  level  by  photon  absorption  exceeds  the  depopula¬ 
tion  by  intermolccul ar  collisions,  second,  diffusion  effects  wherein  the  colli¬ 
sional  effect  is  masked  by  flux  of  ground  state  molecules  into  the  region 
irradiated  by  the  laser.  The  former  of  these  was  treated  in  some  detail  in 
Section  II  for  1  10:1  I^/CO  mixture  at  ten  atmospheres  pressure  irradiated  by  CO 
laser  radiation.  In  that  study,  the  CO  laser  saturation  intensity,  Is,  at 
which  the  CO  absorption  coefficient  decreased  by  30%  due  to  optical  saturation  was 
calculated  to  increase  from  5.6  MW/cm2  at  2000*‘K  to  12  MW/cm2  at  3000*K. 

Although  these  values  are  consistant  with  the  incident  laser  intensities  used 
in  the  flames  experiments,  it  must  be  remembered  that  collisional  deactivation 
for  diatomic  molecules  such  as  CO  occurs  at  a  much  lower  rate  than  for  poly¬ 
atomic  molecules  such  as  l^O,  and  the  Ig  values  calculated  for  CO  are 
probably  low  compared  to  the  actual  value  for  1^0  so  that  the  observed  gap 
between  the  flame  and  plasma  data  for  is  probably  not  ascribable  to 

optical  saturation  involving  collisional  phenomena. 
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The  second  phenomenon  mentioned  was  cited  in  Ref.  42  for  cancelling 
optical  saturation  involving  collisions  by  renewing,  through  diffusion  into  the 
laser  beam,  the  supply  of  ground  state  molecules  at  a  faster  rate  than  through 
collisional  deactivation  of  excited  molecules.  The  beam  diameter  at  which  the 
rate  of  collisional  deactivation  equals  that  of  diffusion  is  given  by  the 
square  root  of  the  ratio  of  twice  the  diffusion  coefficient  to  the  collisional 
deactivation  rate  and  this  diameter  was  calculated  for  H2/H2O  to  range  from 
9.76Um  to  3.951*m  over  the  temperature  range  between  300°K  and  6000*K  at  10  atm 
pressure.  Since  the  smallest  high  power  beam  diameter  in  the  flame  experiments 
was  roughly  80Mm,  it  is  seen  that  diffusion  effects  cannot  influence  the  value 
of  I  in  these  experiments. 

Significance  of  Result  s_ 

The  most  noteworthy  feature  of  the  results  shown  in  Fig.  42  is  their 
magnitude  compared  to  the  calculated  value.  It  is  seen  that  the  former  are 
consistantly  about  an  order  of  magnitude  greater  than  the  latter  throughout  the 
temperature  range  shown.  At  the  same  time,  it  is  seen  that  a  curve  rough  drawn 
through  the  entire  experimental  data  set  is  nearly  parallel  to  the  calculated 
curve.  These  facts  indicate  that  although  the  magnitude  of  the  experimental 
results  is  significantly  higher  than  calculated,  their  temperature  dependence 
is  consistant  with  that  predicted  by  analysis.  Part  of  the  difference  between 
the  analytical  and  experimental  results  may  stem  from  the  fact  that  the  former 
are  derived  from  measurements  averaged  over  23  cm'*  intervals  so  that  signifi¬ 
cant  excursions  from  them  are  expected  for  measurements  made  over  narrower 
intervals  as  in  the  present  case.  The  expected  degree  of  variation  over  a  25 
cm-*  wide  interval  may  be  seen  in  Fig.  8.  The  important  point  to  be  noted 
is  that  the  measured  value  of  ku  is  so  large  as  to  make  HjO  a  very  attractive 
molecular  seed  for  a  laser  energized  rocket  thruster. 


Results  for  D2O 


The  measured  value  of  ku  (T)  is  shown  for  H2/D2O  mixtures  in  Fig.  43. 

As  was  the  case  in  Fig.  42  several  experimental  runs  are  shown.  The  solid  and 
dashed  curves  resulted  from  plasma  and  flame  data  respectively  taken  at  a  D2O 
pressure  of  approximately  0.3  atm  in  13  atm  t^.  The  solid  line-and-circles 
curve  resulted  from  a  flame  measurement  in  which  the  cell  temperature  exceeded 
the  saturator  temperature  so  that  more  confidence  can  be  placed  on  the  DjO 
pressure  given  than  in  the  other  H2/D2O  experiments.  As  in  the  case  of  the 
H2/H2O  experiments,  the  flame  and  plasma  data  at  the  same  water  partial 
pressure  are  seen  to  be  consistant  with  one  another  and  to  provide  data  from 
about  500  to  4500*K.  In  addition,  the  flame  data  taken  at  the  lower  DjO 
partial  pressure  is  quite  similar  in  shape  to  the  other  flame  curve,  but  the 
two  curves  do  not  exhibit  magnitudes  which  scale  as  the  DjO  partial  pressure, 
indicating  that  the  0.3  atm  D2O  partial  pressure  given  is  probably  an  upper 
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bound  to  the  value  actually  present  in  those  experiments.  The  fact  that  the 
maximum  observed  flame  temperatures  do  not  coincide  may  be  due  in  part  to 
experimental  error  at  the  high  temperature  end  of  the  curves  as  discussed  in  a 
preceding  section  and  also  in  part  to  the  fact  that  the  two  sets  of  measurements 
were  taken  in  slightly  different  vertical  positions  in  the  cell  with  the  hotter 
flame  measurement  being  taken  quite  near  the  high  power  laser  beam  focal  spot 
and  the  colder  flame  measurement  being  taken  slightly  above  it. 

There  are  two  analytic  curves  presented  in  Fig.  A3.  The  first  curve, 
labled  "D2O"  was  calculated  for  the  indicated  H2/D20  mixture  in  which  the 
only  effect  of  increasing  temperature  was  to  lower  the  gas  density  and  change 
the  internal  energy  level  population  distribution  in  D2O  concentration  due  to 
the  establishment  of  chemical  equilibrium  which,  for  the  mixture  indicated, 
would  cause  D2O  to  vanish  with  1^0  and  HDO  becoming  the  principle  coupling 
molecule.  These  chemistry  effects  are  included  in  the  curve  labled  "HjO''. 

The  chemical  kinetics  driving  the  H2/D2O  system  to  equilibrium  are  slow  at 
low  temperatures  and  increase  in  speed  with  temperature  so  that  an  experimentally 
measured  ku  (T)  curve  for  H2/D2O  is  expected  to  follow  the  "DjO"  curve  at 
low  temperatures  and  the  "l^O"  curve  at  higher  temperatures  with  the  transi¬ 
tion  between  the  two  curves  occurring  at  the  temperature  at  which  the  time 
needed  to  establish  chemical  equilibrium  equals  the  residence  time  of  the  gas 
in  the  high  temperature  region.  Examination  of  the  dashed  flame  data  curve 
indicates  that  in  this  experiment,  the  noted  condition  occurs  in  the  500-900°K 
interval.  Most  important  however  is  that  the  D2O  data  indicates  this  molecule 
to  be  an  excellent  seed  molecule  coupler  with  the  advantage  over  1^0  of 
having  a  larger  absorption  coefficient  than  1^0  at  lower  temperatures, 
indicating  this  molecule  should  be  easier  to  activate  to  maximum  absorbt i vi ty . 

Results  for  NH^ 


The  ku  (T)  curve  for  the  Hj/NH^  system  is  shown  in  Fig.  AA  with  only 
Che  flame  data  being  shown.  The  ammonia  molecule  is  known  to  couple  very  well 
to  CC>2  laser  radiation  at  room  temperature,  but,  due  to  its  lack  of  chemical 
stability  at  temperatures  above  500*K,  its  usefulness  is  expected  to  be  limited 
to  relatively  low  temperatures.  The  results  shown  in  Fig.  AA  are  therefore 
quite  surprising,  for  instead  of  decreasing  from  its  measured  room  temperature 
value,  shown  as  the  solid  square,  the  value  of  ku  instead  remains  constant 
until  1000'K  is  reached  at  which  point  it  increases  more  than  an  order  of 
magnitude  in  value,  reaching  a  maximum  at  2000°K.  The  reason  for  this  behavior 
is  not  known,  but  one  possibility  is  that  the  spatial  region  over  which  the 
measurements  are  taken  is  small  enough  and  the  chemical  kinetics  slow  enough 
that  the  equilibrium  condition  is  not  reached  until  temperatures  in  excess  of 
2000*K.  In  contrast,  in  the  presence  of  plasma,  the  NH^  concentration  in  the 
cell  is  seen  to  fall  by  about  AO %  in  ten  seconds  when  exposed  to  the  high 
temperatures  associated  with  the  plasma. 
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The  persistance  of  a  large  ku  value  at  high  temperatures  makes  NHj  a  more 
attractive  coupling  candidate  molecule  than  was  heretofore  expected. 

Coupling  Measurements 

Irtt  roduct  ion 

In  order  to  gain  insight  to  the  properties  of  the  laser  sustained  flame  in 
11, /water  mixtures,  measurements  were  carried  out  to  determine  the  degree  of 
couplitg  between  the  high  power  laser  beam  and  these  mixtures  as  a  function  of 
vat' r  concentration  and  beam  power.  In  order  to  do  away  with  the  problem  of 
condensation  within  the  cell,  the  latter  was  kept  at  temperature  higher  than 
the  saturator  temperature  in  each  case.  The  range  of  partial  pressures  was  20 
torr  to  180  torr  for  1^0  and  20  torr  to  130  torr  for  D2O.  In  each  run, 
the  increase  in  the  centerline  temperature  determined  from  the  interferometric 
measurement  was  used  as  a  measure  of  the  coupling.  In  the  l^/l^O  measurements 
the  measurements  indicated  little  if  any  beam-mixture  coupling  over  the  entire 
partial  pressure  range,  and  only  the  results  for  the  H2/D2O  measurements 
will  be  presented.  Of  the  latter,  a  noteworthy  result  was  the  attaining  of  a 
very  hot,  3900'K,  flame  in  only  one  of  the  runs,  at  76  torr  partial  pressure 

and  7. 1  kW  laser  power.  The  ku  (T)  information  obtained  in  this  experiment  is 

shown  in  Fig.  43.  The  observed  temperature  rise  in  the  remaining  H2/D2O 
experiments  in  this  series  was  nearly  a  factor  of  ten  smaller,  even  though  one 
of  these  was  carried  out  at  higher  I^O  partial  pressure,  130  torr,  and  the 
same  beam  power.  However,  another  difference  between  these  two  H2/D2O 
experiments  was  that  the  former  was  carried  out  at  a  cell  temperature  25K" 

cooler  than  the  latter,  so  that  flame  creation,  like  plasma  sustaining,  is 

apparently  critically  dependent  on  temperature  gradients  in  the  high  power  beam 
input  window  and  their  effect  on  the  quality  of  the  beam  to  be  focused  into  the 
gas  mixture. 

The  Measurements 

The  coupling  measurements  consisted  of  noting  the  change  in  the  axial 
temperature  measured  with  the  interferometer.  It  is  re-emphasized  that  in  only 
one  of  the  experiments  reported  on  in  this  section  was  a  laser  sustained  flame 
attained  presumedly  because  of  the  input  window  problem.  The  remaining  measure¬ 
ments  are  best  understood  by  modeling  the  laser-beam-mixture  interaction 
region  as  a  heat  source  of  dimension  R  cooled  by  bouyancy  driven  gas  flow.  The 
balance  between  the  absorption  of  energy  from  the  incident  laser  beam  of 
intensity  I  by  the  gas  with  absorption  per  cm  a  resulting  in  a  center  line 
temperature  rise  AT  cooled  over  the  dimension  R  by  gas  flow  of  speed  V  is  given 
by 

PC  pVAT 

Ia  R  ( 3U> 
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If  I  is  given  by  the  ratio  of  beam  power  p  to  its  area  A,  and  a  is  approximated 
for  D2O  in  the  temperature  range  from  300  to  1000°K  by 


a  =  a0+  AT  da 


(35) 


Then  the  quantity  a^p/AT  is  a  linear  function  of  P  the  incident  beam  power 
and  can  be  written  in  the  form 


^-=a~'-P(Ja  (36) 

where 

a  =  -  R-  - 

AV^jCp 

The  absorption  per  cm  of  the  mixture  was  measured  under  pre-run  conditions  to 
be  0.026  cm  *  at  130  torr  D2O  pressure  and  368°K,  and  this  value  was 
combined  with  the  experimental  values  of  P  and  AT  to  calculate  aQ P/AT  as  a 
function  of  P  as  shown  in  Fig.  45.  The  results  denoted  by  the  solid  points  are 
seen  to  be  consistant  with  this  simple  model,  giving  values  of  1.37  cm  deg/W  and 
5.3  x  10  ^  cm  *  deg-*,  for  a  and  5a  respectively.  The  lone  flame  results, 
denoted  by  the  cross,  is  seen  to  fall  far  from  the  others  due  to  the  much  larger 
AT  value  attained  in  this  experiment. 

In  studying  the  effect  of  D2O  concentration  on  the  coupling  Eq.  (3*0  was 

solved  for  using  the  value  of  a  and  derived  from  Fig.  45,  and  the 

value  of  a  so  derived  from  the  data  is  plotted  vs.  the  D~0  concentration 
o  r  2 

in  Fig.  46.  Within  rather  large  error  limits  it  is  seen  that  a  straight  line 
which  passes  through  the  origin  can  be  drawn  through  the  data.  The  slope  of 
this  line,  0.23  cm-*  amagat-  is  the  experimentally  derived  value  of  k  for 
I^O  in  the  temperature  region  near  330°K  and  is  significantly  larger  than  the 
calculated  value  of  C  plotted  for  D2O  in  Fig.  26. 
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SECTION  V 
DISCUSSION 


Summary  of  Results 

An  analytical  and  experimental  program  has  been  carried  out  to  evaluate 
the  usefulness  of  binary  mixtures  of  hydrogen  with  H20,  D20  and  NHj  for 
coupling  energy  transported  from  a  distant  source  by  a  laser  beam  to  heat  a 
rocket  working  fluid.  Specifically,  the  desired  information  was  the  tempera¬ 
ture  dependence  of  the  absorption  per  cm  of  these  mixtures  for  temperatures  as 
high  as  6000*K.  To  obtain  this  information,  the  gas  mixture  under  investiga¬ 
tion  was  placed  in  a  high  pressure  cell  in  the  presence  of  a  hot,  1 0000-1 5000*K, 
CW  laser  sustained  plasma  with  the  cell  in  turn  placed  in  a  diagnostic  apparatus 
which  measured  simultaneously  the  absorption  per  cm  and  the  temperature  at 
points  about  the  axis  of  the  laser  beam  sustaining  the  plasma.  These  measure¬ 
ments  indicated  that  the  absorption  per  cm  for  the  H2/H20  and  H2/D20 
mixtures  significantly  exceeds  the  magnitude  expected  from  the  analytic  study 
while  also  exhibiting  a  temperature  dependence  consistant  with  the  establishment 
of  chemical  equilibrium  within  the  region  of  elevated  temperature.  The  magni¬ 
tude,  0.1  to  1.0  cm-*  in  the  1000-5000*K  temperature  range,  makes  either  of 
these  mixtures  extremely  attractive  for  use  in  a  laser  energized  rocket  thruster. 

In  addition,  in  each  of  the  mixtures  studied  it  was  found  that  the 
interaction  between  the  focused  high  power  laser  beam  and  the  mixture  was 
often  intense  enough  to  create  a  hot  luminous  region  at  the  focal  spot, 
referred  to  in  this  report  as  a  laser  sustained  flame.  Although  the  center 
line  temperatures  in  the  flame  were  quite  high,  as  high  as  3900*K  in  one  case, 
the  flames  were  cooler  than  the  plasma  and  served  to  supplement  the  plasma  data 
by  extending  the  observed  temperature  range  to  lower  values.  The  magnitude  of 
the  absorption  per  cm  measured  in  the  l^/water  flames  was  consistant  with 
chat  in  the  plasma  experiments  even  though  the  high  power  laser  intensity  in 
the  measurement  region  was  much  higher,  indicating  that  optical  saturation 
effects  are  not  important  for  the  H2/water  mixtures  and  CW  C02  laser 
intensities  up  to  10  MW/cm^  used  in  these  experiments.  For  the  H2/NH-j 
flame,  the  measurements  indicated  the  NH-j  molecule  to  be  absorbing  at  tempera¬ 
tures  as  high  as  3000*K  in  contrast  to  the  prediction  of  the  destruction  of 
NH^  at  chemical  equilibrium  for  temperatures  above/500*K.  This  indicates 
that  the  kinetics  driving  the  Hj/NH^  system  to  equilibrium  may  be  slow 
enough  to  permit  use  of  ammonia  in  a  rocket  thruster  at  temperatures  higher 
than  previously  believed  possible.  In  the  plasma  experiments  involving  ammonia, 
the  destruction  of  this  molecule  at  the  temperatures  associated  with  the  plasma 
was  very  apparent,  indicating  that  there  is  a  limit  imposed  by  chemical  effects 
on  the  usefulness  of  NH^. 
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Unknown  Flame  Properties 

One  of  the  more  interesting  results  of  this  study  is  the  observation  of 
the  ability  of  the  high  power  laser  beam  to  raise  the  temperature  of  the 
mixture  under  study  to  a  level  useful  in  a  laser  energized  rocket  thruster. 

Very  little  is  known  about  the  ease  with  which  such  flames  are  made.  It  is 
known  that  for  a  focal  spot  diameter  of  about  0.01  cm,  flames  can  be  created  in 
H^/DjO  and  l^/NH^  mixtures  with  7  kW  of  laser  power.  However,  there  is 
evidence,  described  in  Section  IV,  that  if  the  laser  beam  quality  is  decreased, 
e.g.,  by  distortion  on  passing  through  the  cell  window,  and  the  focal  spot 
diameter  becomes  larger,  creation  of  the  flame  is  more  difficult.  This  problem 
was  also  encountered  in  sustaining  plasma  with  laser  radiation,  and  one  is 
therefore  led  to  the  conclusion  that  just  as  the  laser  power  needed  to  create 
and  sustain  plasma  in  gases  exhibits  a  beam  diameter  dependence,  J  so  will  the 
power  needed  to  create  and  sustain  a  flame.  Put  another  way,  it  is  not  clear 
that  just  because  7  kW  focused  into  a  spot  10  cm  in  area  will  create  a 

•  9  t 

flame,  7  MW  focused  into  a  spot  1  cm  in  area  will  do  the  same  even  though 
the  intensity  is  the  same  in  both  cases.  Other  questions  of  interest  concern 
the  beam  diameter  at  which  the  time  needed  for  NH^  destruction  is  exceeded  by 
the  diffusion  lifetime  of  NH^  inside  the  beam,  and  the  possibility  that  the 
beam  power  needed  to  create  a  flame  in  H2/H2O  might  be  lowered  by  the 
addition  of  a  NHj  "Priming  charge"  to  bring  the  HjO  absorption  per  cm  high 
enough  to  permit  ignition  of  the  ^0  flame. 

In  any  event,  the  discovery  of  the  H2O  laser  sustained  flame  was  an 
unexpected  result  of  this  Air  Force  sponsored  research  program  and  one  which 
might  prove  extremely  valuable  in  the  area  of  laser  energized  rocket  thrusters. 
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CARBON  MONOXIDE  LASER  LIGHT  SOURCE 
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Tin*  carbon  monoxide  laser,  wliijh  was  to  he  used  a;;  the  diagnos!  i(  1  ij'lit 
source  at  fre<|uencies  near  2000  cm  ,  was  built  largely  based  on  tin-  worl 
described  in  Ref.  44.  The  laser  cavity  it  sell  was  approximately  2  ’  met.  is 
in  length,  with  a  10  m  radius  of  curvature  gold-coated  copper  high  reflect  ivif 
mirror,  and  a  flat  germanium  95  percent  reflectance  mirror  output  coupler.  The 
active  medium  was  an  electric  discharge,  powered  by  current  regulated  power 
supplies,  in  a  gas  mixture  of  approximately  16  torr  pressure  containing  hilium, 
carbon  monoxide,  nitrogen,  and  oxygen  in  the  ratio  of  1500  :  40  :  40  :  1  .  The 
helium  was  passed  through  a  coil  of  copper  tubing  immersed  in  liquid  nitrogen 
to  freeze  out  any  water  in  the  helium  prior  to  entering  the  discharge  region. 

The  discharge  tube  had  two  identical  sections  made  from  py rex  and  joined  in 
the  middle  by  a  short  pvrex  cross  which  held  the  gas  outlet  as  well  as  the 
common  cathode  to  each  section’s  discharge.  The  inside  diameter  of  the  discharge 
tubes  was  3/4  in.  The  active  medium  was  sealed  at  each  end  with  water  cooled 
calcium  fluoride  windows  mounted  at  Brewster's  Angle,  and  the  region  between 
each  window  and  cavity  mirror  was  continually  flushed  with  dry  nitrogen  to 
drive  out  water  vapor  which,  if  left  in  the  cavity,  absorbed  cavity  radiation, 
thereby  becoming  thermally  and  optically  inhomogeneous  and  causing  the  laser 
output  to  be  both  temporally  unstable  and  diminished  with  respect  to  both  output 
and  beam  quality.  The  laser  could  be  operated  both  at  ambient  temperature, 
with  water  cooling,  and  cryogenic  temperature  with  liquid  nitrogen  cooling. 

The  cooling  jackets  for  each  discharge  section  were  approximately  75  cm  long 
with  a  metal  bellows  blown  into  the  middle  of  the  outer  wall  of  the  jacket. 

The  purpose  of  the  bellows  was  to  relieve  stress  caused  by  the  difference  of 
length  which  was  assumed  by  the  inner  wall  of  the  cooling  jacket  at  liquid 
nitrogen  temperature  compared  to  the  outer  wall  which,  of  course,  assumed  some 
temperature  between  liquid  nitrogen  and  ambient.  Experience  showed  that  cracking 
of  the  jackets  routinely  occurred  in  the  absence  of  the  bellows.  The  inter-wall 
distance  in  the  jackets  was  1.6  cm.  It  was  found  that  such  a  large  distance  was 
necessary  to  have  the  discharge  tube  completely  immersed  in  the  liquid  nitrogen 
while  leaving  enough  free  space  at  the  top  of  the  jacket  to  allow  nitrogen  vapor 
to  collect  and  exhaust  without  causing  wave  formation  in  the  liquid  nitrogen  which 
resulted,  at  smaller  Inter-wall  spacing?,  in  blowing  liquid  nitrogen  out  of  the 
exhaust  port.  With  water  cooling,  the  multiline  output  of  the  laser  was  as  high 
as  20  watts,  with  13  watts  being  obtained  routinely  over  hour-long  periods.  With 
liquid  nitrogen,  70  watts  of  multiline  output  power  were  obtained. 
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"STOR2" 


ft 


2  REM  N  IS  ARRAY  TRACE  NUMBER — NUMERIC 

3  LET  N-e 


4  REM  R#  IS  RUN  NUMBER— SPELLED  OUT 

5  PRINT  " INPUT  RUN  NUMBER-  THEN  •RETURN’" 

MPUT  R* 

LQAO  "DPO" 

ATTACH  §2  AS  DPO 
WHEN  #2  HAS  "CBl"  COSUB  46 
RETURN 

PUT  "PI"  INTO  42 
PACE 

LET  N*-STR<N> 

LET  RN*-R*4Nf 

OPEN  #1  AS  DXJ  RNf  FOR  WRITE 

PRINT  INPUT  WHUEKJRM  T'J  UEkWLAL  rLUC-iN-  HP  TEFQRT1S  MUST  SPAN  10  OIUS 
PRINT  "UNGROUND  PROBE  AND  ADJUST  KNOBS  FOR  PROPER  SENSITIVITY" 

PRINT  “PRESS  'RETURN' WHEN  READY 
WAIT 
PAGE 

WALEFORM  WA  IS  AC311  LW-HA4-UA4 
PUT  "STO"  INTO  42- "A" 

PUT  "MEM"  INTO  42 
PUT  "ZHO"  INTO  42- "A" 

PRINT  "TO  ACQUIRE  GROUND-REFERENCE  LEVEL-  GROUND  PLUG-IN" 

"PRESS  PROG  ChLL*  BUTTON  14  ON  DPO-PRESS  'RETURN'  WEN  READY" 


PRINT 

WAIT 

PAGE 

PRINT 

PRINT 

WAIT 


"TO  ACQUIRE  AUE  MALCFQRM,  UNGROUND  PLUG-IN" 
"PRESS  RETURN'  WEN  READY- 


PUT  "STO"  INTO  42- "A" 

PUT  "MEN"  INTO  42 

GET  HA  FROM  42- "AUE20- A" 


10 r) 


6  OIH  MEK51 1  > 

9  OIH  MA(511),MZ<911) 

10  INPUT  Zt 

20  OPEN  •!  AS  0X1  Z*  FOR  READ 
90  READ  tl,UZ 
40  CLOSE  *1 
90  INPUT  A* 

60  OPEN  tl  A3  0X1  A*  FOR  READ 

70  READ  01,  HA 

71  CLOSE  tl 


100  > 


90  GRAPH  MZ<0> 

100  WAIT 
110  PAGE 

120  GRAPH  MAC  0 1 100) 
130  MAIT 
140  PAGE 

190  I»*UT  A/B'C/O/E 
160  FOR  I -A  TO  B  STEP 
170  MKI>-UA(I*0>»€ 
180  t€XT  I 
109  GOTO  3900 
190  PAGE 

200  GRAPH  MB, HZ 

201  MAIT 


210  IfWT  J 

IF  >0  THEN  GOSUB  1000 


230  OPEN  tl  AS  OKI  B*  FOR  MRITE 
240  MRITE  01  ,m 
CLOSE  01 


INPUT  OE 


108 


1013  PAGE 
1020  GOTO  280 
1030  RETURN 


2010  GRAPH  HB 

2020  INPUT  Bt 

2030  RETURN 

3000  GRAPH  MB< 360  309  > 

3010  WAIT 


INPUT  IB'UB 
FOR  I«0  TO  311 

leu  >-»a  >-u8*coa  tti 

NEXT  I 
GOTO  190 


14139*<IB-I>/167  3>> 


READY 


* 


le  DIM  A<311  >,B<511  >,C<311  ),D<51 
20  INPUT  At,B#,C#,Df 
30  OPEN  tl  AS  DX1  At  FOP  READ 
40  READ  il » A 
30  CLOSE  tl 

60  OPEN  tl  AS  0X1  B#  FOR  READ 
70  READ  tl  >  B 
00  CLOSE  tl 

90  OPEN  tl  AS  0X1  Ct  FOR  REPO 

100  READ  tl,C 

110  CLOSE  tl 

120  PACE 

130  GRAPH  A;B/C 

139  WAIT 

140  0-<  A-B-C  V<  2*SQR<  B*C  )  ) 

130  PAGE 

160  GRAPH  0 

161  WAIT 

162  PAGE 
170  INPUT  M 

100  IF  M-0  THEN  GOTO  300 

190  OPEN  tl  AS  0X1  Df  FOR  l*ITE 

200  WRITE  tl,0 

210  CLOSE  tl 

300  END 


READY 

* 
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le 

13 

16 

20 

90 

40 

43 

30 

60 

63 

66 

70 

00 

89 

90 

91 


DIM  A<  300  )/B<  306  ).MP<31 1  ),PX<  311) 
PI-3  14139 
X— PI 

FOR  1-0  TO  500 
<tt  I  >-X 
0<  I>»COS<X> 

X-X*PI^300 
NEXT  I 
INPUT  R* 


AX-0 
OPEN  #1 
REPO  il 
PAGE 
CLOSE  41 
GRAPH  HP 


PS 

UP 


0X1  R*  FOR  REPO 


92  HP IT 

100  INPUT  MN.-MX 

120  FOR  I"4«4  TO  MX 

121  IF  P8SC  MAC  I  >  X 1  THEN  GOTO  127 

122  IF  MAC  I  »0  THEN  GOTO  123 

123  AXC  I  >— PI 

124  GOTO  170 
123  PKC  I  >«0 

126  GOTO  170 

127  FOR  J-0  TO  300 

130  IF  P0S<  HP<  I  )-B<  J )  X-3E-03  T1CN  GOTO  136 
140  NEXT  J 
130  PX<  I  >-A<  J) 

170  NEXT  I 
170  PAGE 
179  PRINT  Rl 

100  GRAPH  UP, AX 

101  MAIT 


r  * 


190  INPUT  St 

200  OPEN  il  MS  DX1  St  FOR  WRITE 
210  WRITE  il /AX 
220  CLOSE  »1 
230  EM) 


READY 

* 


1  REM  OATH  SMOOTHING  ROUTINE  "HASAG" 

9  DIM  MAC  51 1  ),HB(511  )-MC(511  )»MD(311  )/A(9/ 10) 

10  IfFHJT  R#/$$ 

20  GOSUB  2000 

21  GRAPH  MB 

22  WAIT 

24  GRAPH  HA 

25  INPUT  MN/MX 

90  GRAPH  HAcNN  NX) 

40  INPUT  SA 

50  IF  SA-1  THEN  110 

51  MC«0 

52  M>>0 

€0  INPUT  JA 
70  GOSlff*  3008 
80  OISPLAY  1  > MCC HN : MX  > 

90  INPUT  SA 
100  GOTO  50 
110  MA-0 
in  w-e 
113  GOSUB  4000 

2000  OPEN  tl  AS  0X1  Rf  FOR  READ 

2001  EOF  *1  GOTO  2020 
2802  READ  41 /HA 
2003  CLOSE  01 

2025  GOSUB  5000 
2030  RETURN 

3000  FOR  J-NN  TO  MX  STEP  JAM 

3001  MCC  J>-0 

3002  MDC  J>0 


FOR  I-J  TO  J*JA 
O+CUH** I ) 

NEXT 


2 


3809  MCC  J+JfV2  >-04":  JHM  ) 

3810  HOC  J-*- J4V2  >»WB< J+ JA/2 ) 

3820  NEXT  J 

3890  RETURN 
4800  REM 
4820  I 1*0 

4090  FOR  I-MN  TO  MX 
4040  IF  WC<I»e  THEN  4000 
4090  MA<  1 1  >-WC<  I  ) 

4060  HBC  1 1  >-WCK  I  > 

4061  PRINT  UMlI>.te(n>,U 
4070  II-IH-1 

4000  rcxT  I 
4090  INPUT  K0>  JJ 

4891  REM  ZtJJ  IS  THE  LARGEST  POtCR  OF  Y 

4092  REM  JJ<-6  AS  DEFINED  BY  DIMCA) 

4093  A*0 

4108  FOR  l-l  TO  JJ*t 
4110  FOR  J-l  TO  JJ+2 

4120  FOR  K-K0  TO  II-l 

4121  IF  J-J-H2  THEN  4141 

4190  A<  I » J  )-A<  I  •  J  HABS<  HB(  K  )  X  2t<  1  ♦  J-2  >  ) 
4132  GOTO  4190 

4141  A<  I,  J>-A<  1/ J>H#XK)*AeS<l<KK>  W  2*1-2  > 
4198  NEXT  K 
4168  »CXT  J 

4178  *CXT  I 

4179  Ml-JJ+l 

4179  FOR  1-1  TO  JJ+1 

4180  >tt-ABS<A(I,J)> 

4280  L-I 

4210  KK-IM 

4220  FOR  K-KK  TO  Ml 

4290  IF  XN>A8S<AOC,I>>  THEN  4260 

4248  )!Wi8((XKll)) 


122 


4256  L«* 

4260  hCXT  K 

4270  FOR  J=1  TO  h*Kl 

4200  B-AC  I ,  J  ) 

4290  <*  I,  J>-ACL,  J) 

4300  ACL,  J>-8 
4310  NEXT  J 
4320  H-NN+1 

4330  Aa.MHKI.MVfXLl) 

4340  IF  Hc*I  THEN  4370 

4350  M-M-l 

4360  GOTO  4330 

4370  IF  I<-1  THEN  4400 

4300  hHl 

4390  GOTO  4440 

4400  N-l 

4410  N-N+l 

4420  IF  N»I  T»CN  4410 
4438  IF  N>NN  THEN  4490 

4450  I>«*I«H> 

4460  IF  HOI  THEN  4410 
4470  rWT-1 
4460  GOTO  4450 
4490  »CXT  I 
4500  PRINT  A< 

4509  INPUT  MK 

4510  FOR  I-4C0  TO  HK 
4520  MCC  I  »e 

4530  FOR  J»1  TO  JJ+\ 

4540  MC<  I  >^C<  I  >*0<  J,JJ+2 >tP0S< X  )  W 2t>2 > 

4550  NEXT  J 

4551  PRINT  HC<  I  >/MBC  I  > 

4560  NEXT  I 

4565  PRINT  R#i-  US  "i80>"  NN»")»0«i- 


4579  WAIT 

4575  GRhPH  WC(MN  MX) 

4569  WAIT 

4563  GOSUB  5666 

4585  FOR  I-t*  TO  MX 

4596  MC<  I  >-6 

4595  FOR  J-l  TO  JJM 

4686  MC<  I  >-MC<  I  J,  J  J*2  >*ABS<  HBU  >  W  2tJ-2  > 

4665  NEXT  J 
4610  NEXT  I 

4624  PRINT  "MAKE  SURE  RIGHT  OISC  IS  IN  0X1 M 

4625  WAIT 

4626  GOSUB  2666 

4627  PRINT  A< 1  * JJ*1 , JJ*2 > 

4628  MC<213>-AU, JJ«-2> 

4636  GRAPH  M*nN  =  HX>,MC<»t4  NX) 

4635  WAIT 
4646  END 

5666  FOR  1-511  TO  6  STEP  -1 
5616  WUj-5  35E-63*<I-213> 

5626  NEXT  I 
5636  RETURN 
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3  dim  r<  iei  >, «l<  iei  >,  n  tei  >, tt<  lei 

6  OIM  N(7>,P<7>,C<7,7>,L<7>,M<7> 

7  OIM  A0<7> 

10  INPUT  T0,P8,PP,B,RA,P,RP 

11  B0-8 

13  INPUT  RN.RX 

17  PI-3  14159 

18  K8>-*N 
FOR  1-8  TO  18 
H-8 
Z-RB 

X-SOK  Z-2-R<  I  V2> 

COSUB  1888 
FOB  J-l  TO  7 
L<  J>— X^<2*J-1  VP  I 

NEXT  J 
GOSUB  2888 
Z"ffP 

X-SQRC  ZtZ-fK  I  V2  > 

FOB  J-l  TO  7 
HC  J>— X~<2*J-1  VP  I 

NEXT  J 
GOSUB  1888 
GOSUB  3888 

T<  I  )-2  68E-03JN+3  8E-83*P0/T0 
GOSUB  4888 

R<  !♦!  I  >K RX-RN  V10 

A»BA 
NEXT  I 

PBINT  T8,Pe,PP,A,RA,P,RP 
FOB  1-8  TO  18 

PBIMT  RU>,RL<I),T<I>,TT<I> 
NEXT  I 
MBIT 


12b 


r  f 


63  PAGE 

65  GRAPH  AL<  0  10  > 

66  WAIT 

67  PAGE 

70  GRAPH  1(0  10  > 

72  WAIT 

73  PAGE 

75  GRAPH  TT<0  10) 

80  WAIT 
85  END 

1000  C-0 

1001  IF  H-0  THEN  1008 

1002  A-P 

1006  n>i 

1003  FOR  J-l  TO  7 

1010  C<  J » 1  >-Z~<  2*  J-2  >*0(  J  >*2*J 
1015  IF  J-l  THEN  GOTO  1612 
1020  IF  J-2  THEN  GOTO  1522 
1030  IF  J-3  THEN  GOTO  1432 
1040  IF  J-4  THEN  GOTO  1342 
1090  IF  J-5  THEN  GOTO  1292 
1060  IF  J-6  THEN  GOTO  1162 

1072  C<7,2>— 4*2-10*0(7  >*14 

1073  C(7,3*8*Z-8*0(7)*14 

1074  C<  7  #  4  >— 64'7*2-6*0(  7  >*14 

1075  C<7,5>-384/63*Z-4*0<7>*14 

1076  C<7,6>— 512/231*Z-2*0(7>*14 

1077  C<  7. 7  >-3072/  9009 *0(  7  X14 
1076  GOTO  1612 

1162  C(  6, 2  >— 10^3*2 -8*0(6  >*12 

1163  C(6,3)-16x3*2-6*0(6>tl2 

1164  0(6, 4  >-32/7*2-4*0(6  >*12 
1169  C(  6, 5  >-128/63*2-2*0(6  >*12 

1166  C< 6, 6 >—256/693*0(6 X12 

1167  GOTO  1612 


127 


1232  C<3,2>— 8^3»2-6*A<5)fie 

1233  C<3,3>-16/'5*2M:WK5>*10 

1234  C<5,4>— «4'33*Z-2ttK3>tl0 

1235  C(3/5>-128/'313*a(3>*t0 

1236  GOTO  1612 

1342  CC4,2>— 6/3»2M*A<4)*8 

1343  C<  4 , 3  >-8/3*Z-2*A<  4 >*8 

1344  tt4*4>a-16'33*K4)«0 

1345  GOTO  1612 

1432  C<  3 , 2  >»- 4/3*2~2*X  3  >*6 

1433  C<3,3)-8/15*A<3>*6 

1434  GOTO  1612 

1322  C<2,2>— 2/3*A(2>*4 

1323  GOTO  1612 

1612  NEXT  J 

1613  RETURN 
2008  *L<I>-e 
2001  FOR  J-l  TO  7 
2010  FOR  K-l  TO  7 

2020  RL<I>-AL<l>H.<K)|tC<XK> 

2030  NEXT  K 
2040  1CXT  J 
2050  RETURN 

3000  H*0 

3001  FOR  J-l  TO  7 
3010  FOR  K-l  TO  7 
1020  N-N*fKK)*C(  J,K> 

3030  NEXT  K 
3040  1CXT  J 
3090  RETURN 

4000  TT<  I  >-6 . 277E-03*PP/T<  1 ) 

4010  IF  TT<  1X3000  THEN  GOTO  4090 
4020  TT< I )-4  965E-03*FP'T< I ) 

4090  IF  TT<  I  >>4500  THEN  GOTO  4090 
4040  TT<  I  )-0 


120 


